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PREDICTION OF SUBSONIC VORTEX SHEDDING 
FROM FOREBODIES WITH CHINES 

Michael R. Mendenhall 
Daniel J. Lesieutre 

Nielsen Engineering & Research, Inc. 

SUMMARY 

An engineering prediction method and associated computer code VTXCHN to predict 
nose vortex shedding from circular and noncircular forebodies with sharp chine edges in 
subsonic flow at angles of attack and roll are presented. Axisymmetric bodies are 
represented by point sources and doublets, and noncircular cross sections are transformed 
to a circle by either analytical or numerical conformal transformations. The lee side vortex 
wake is modeled by discrete vortices in crossflow planes along the body; thus the three- 
dimensional steady flow problem is reduced to a two-dimensional, unsteady, separated 
flow problem for solution. Comparison of measured and predicted surface pressure 
distributions, flow field surveys, and aerodynamic characteristics are presented for 
noncircular bodies alone and forebodies with sharp chines. 


INTRODUCTION 


Current flight vehicle applications requiring high aerodynamic performance can 
involve a variety of noncircular body shapes in subsonic flow at high angles of attack and 
nonzero roll angles. When these bodies have sharp edges, chines, or wing leading edge 
extensions, for example, separation is fixed at the sharp edge, and the lee side vortex wake 
is different from the traditional wake formed by smooth body separation vortices. As in the 
case of smooth bodies, the chine-body vortex shedding characteristics are directly 
influenced by the body cross-sectional shape and the local flow conditions. It is desirable to 
model the lee side vortex wake by means of a rational method capable of considering a 
variety of body shapes over a wide range of incidence angles and Mach numbers up to the 
critical speed. It is important that the separation vortex wake-induced effects on the 
nonlinear aerodynamic characteristics of the configuration be handled properly with a 
method which correctly represents the physical characteristics of the flow field 
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The phenomena of interest are the sheets of vorticity formed when the fluid flow 
separates from the sharp edges on both sides of the body (Fig. 1). At moderate angles of 
incidence, the vorticity rolls up into a symmetrical vortex pair, but at higher angles, the 
vorticity becomes asymmetric and bursting occurs. These latter phenomena are beyond the 
scope of the current method presented herein. 

A method to predict vortex shedding flow phenomena from smooth bodies, both 
circular and noncircular bodies in subsonic flow, is described in References 1 and 2. The 
extension of the prediction method to supersonic flow is described in Reference 3. The 
purpose of this report is to describe an engineering prediction method and associated 
computer code developed to calculate the nonlinear aerodynamic characteristics and flow 
fields of noncircular bodies with sharp edges at moderate angles of attack at speeds up to 
the critical speed. This prediction procedure is based on an extension of the previous 
methods for smooth bodies (Ref. 1). 

The objectives of the method are to use a three-dimensional, attached flow, potential 
method to represent the body and a two-dimensional, incompressible, separated flow model 
to calculate the lee side vortex shedding from the body alone at angle of attack and angle of 
roll. The predicted pressure distribution on the body under the influence of the free stream 
and the separation vortex wake is used to calculate the aerodynamic loads on the body. 

Conformal mapping techniques are used to transform noncircular cross sections to a circle 
tor calculation purposes. 

The following sections of this report include a discussion of the approach to the 
problem and a description of the analysis and flow models required to carry out the 
calculation. The prediction method is evaluated through comparison of measured and 
predicted results on a variety of body shapes, including elliptical cross sections and 
forebodies with chines. A user's manual for the computer code is also included. The 
manual consists of descriptions of input and output and sample cases. 
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LIST OF SYMBOLS 


a horizontal half axis length of elliptic cross section 

A k coefficients of conformal transformation, Eq. (8) 
b vertical half axis length of elliptic cross section 

c f skin-friction coefficient, Eq. (71) 

c p normal-force coefficient per unit length, Eq. (61 ) 

c y side-force coefficient per unit length, Eq. (65) 

C A axial-force coefficient, Eqs. (70) and (74) 

Cj rolling-moment coefficient, Eq. (69) 

C m pitching-moment coefficient, Eq. (63) 

C n yawing-moment coefficient, Eq. (67) 

Cp pressure coefficient, Eq. (53) 

C Y side-force coefficient, Eq. (66) 

C N normal-force coefficient, Eq. (62) 

D diameter of circular body, or diameter of noncircular body with equivalent cross 

section area, D = 2r eq 

G complex velocity component, Eq. (38) 

K total number of Fourier coefficients used to describe transformation, Eq. (8) 

1 f reference length 

L model length 

M x rolling moment about the x-axis 

M pitching moment about the v-axis 

M yawing moment about the z-axis 

M free-stream Mach number 

00 

N normal force 
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p local static pressure 

p OT free-stream static pressure 

free-stream dynamic pressure, (1 /2) p 
Q source strength 

r radial distance from a vortex to a field point 

r' radial distance to a point on a noncircular body, Fig. 4 

r c vortex core radius, Eq. (76) 

r eci equivalent radius of cross section, (S c / 7 r) 1/2 

r Q radius of circle in transform plane 

Re^ Reynolds number based on boundary layer run length and minimum pressure 
conditions, U m £/v 

S reference area 

S c cross sectional area 

u,v,w velocity components in real plane 

u r axial perturbation velocity component from body, Eq. (10) 

U local velocity 

v R radial perturbation velocity from body, Eq. (11) 

Vg vortex-induced velocity, Eq. (75) 

free-stream velocity 
W complex potential, Eq. (26) 

x,y,z body coordinate system with origin at the nose: x positive aft along the model axis, 
y positive to starboard, and z positive up 

x m axial location of center of moments 

Y side force 

a angle of attack 
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a. angle between free-stream velocity vector and body axis 

0 angle of sideslip; also polar angle in a-plane, Fig. 4; also (1-M^ 2 ) 1 /2 

0' local slope of body surface, Fig. 4 

7 ratio of specific heats 

7 r exterior angles of body segment for numerical mapping, Eq. (6) 

T vortex strength 

Ax axial length increment 

5 complex coordinate in an intermediate plane. Fig 2(b) 

t? vertical coordinate in an intermediate plane. Fig. 2(b) 

9 polar angle in v-plane. Fig. 2(a) 

v complex coordinate in circle plane, Fig. 2(a); also kinematic viscosity 

1 run length, Eqs. (58) and (59); or lateral coordinate in an intermediate plane. Fig. 
2(b) 

p free-stream density 

a complex coordinate in real plane, Fig. 2(a) 

r,X lateral and vertical coordinates in circle plane. Fig. 2(a) 

<p roll angle 

<$ velocity potential in real plane 

\p stream function in real plane 

Subscripts and Superscripts 
(“) conjugate of complex quantity 

( ) vortex m 

( ) center of pressure 

( Y incompressible quantity; or surface values in Fig. 4 
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GENERAL APPROACH 


Smooth forebodies of missiles and aircraft at high angles of attack exhibit distributed 
v orticity fields on their lee side due to boundary layer fluid leaving the fuselage surface at 
separation lines. One approach for modeling these distributed vorticity fields has involved 
the use of clouds of discrete potential vortices. Underlying the basic approach is the 
analogy between two-dimensional unsteady flow past a body and the steady three- 
dimensional flow past an inclined body. The three-dimensional steady flow problem is 
reduced to the two-dimensional unsteady separated flow problem for solution. Linear 
theory for the attached flow model and slender body theory to represent the interactions of 
the vortices are combined to produce a nonlinear prediction method. The details of the 
application of this approach to the prediction of subsonic flow about smooth circular and 
noncircular forebodies are presented in References 1 and 2. Other investigators have used 
this approach to successfully model the subsonic flow phenomena in the vicinity of circular 
cross section bodies (Ref. 4). In spite of the pessimistic outlook, Reference 5 provides a 
comprehensive review of many of the discrete vortex methods currently available. 

The purpose of this report is to document the extension of the subsonic analysis of 
Reference 1 to predict the vortex shedding characteristics of forebodies with sharp comers 
or chines. The code VTXCHN was assembled to accomplish the calculation of subsonic 
vortex shedding from these forebody shapes. 


The calculation procedure for smooth bodies is described in detail in Reference 1. The 
procedure for chined forebodies is carried out in a similar manner. The volume of the body 
is represented by discrete point sources and doublets, and the strength of the individual 
singularities is determined to satisfy a flow tangency condition on the body in a 
nonseparated uniform flow at angle of incidence and roll. Compressibility effects on the 
body are included by a Gothert transformation which keeps the cross section shape 
unchanged but stretches the axial body coordinate. Starting at a crossflow plane near the 
body nose, the pressure distribution on the body is computed using the full compressible 
Bernoulli equation. The bpdy shape determines the location of separation as the assumption 
is made that separation occurs at the sharp edge, and the vortex sheet originates at that 
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point. At the separation points, incompressible vortices with their strengths and positions 
determined by the requirement that the Kutta condition at the sharp edge be satisfied are 
shed into the flow field. The trajectories of these free vortices between this crossflow plane 
and the next plane downstream are calculated by integration of the equations of motion of 
each vortex, including the influence of the free stream, the body, and other vortices. The 
pressure and trajectory calculations are carried out by mapping the noncircular cross section 
shape to a circle using either analytical or numerical conformal transformations. The 
vortex-induced velocity contribution to the body tangency boundary condition includes 
image vortices in the circle plane. 

At the next downstream crossflow plane, new vortices are shed, adding to the vortex 
feeding sheet and cloud representing the wake on the lee side of the body. This procedure is 
carried out in a stepwise fashion over the length of the forebody. The details of the 
individual methods combined into the prediction method for chines are described in the 
following sections. 


METHODS OF ANALYSIS 
Background 

Prediction of vortex shedding from configurations with sharp edges at high angles of 
attack has been an important research area for a number of years because of the dominance 
of the vortex field on the nonlinear aerodynamics of fighter aircraft and missiles. Most of 
this work emphasized sharp edged delta wings because of the availability of experimental 
data and the increasing understanding of the flow phenomena. Unfortunately, forebodies 
with chines have not experienced this same intensity of study; therefore, the base of 
knowledge of the flow phenomena is not as well developed. 

As part of the current effort, a brief examination of the methods and techniques applied 
to the prediction of vortex shedding from delta wings and wing-bodies at high angles of 
attack was made to help put the chine problem into perspective. It is not the purpose of this 
work to review all the theoretical analyses of delta wings, others have already accomplished 
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a detailed review (Refs. 5 and 6). It is the purpose of this work to develop an engineering 
prediction method which does not have some of the limitations and restrictions inherent in 
some of these available methods, but one which still provides a preliminary design and 
analysis capability with reasonable accuracy and economy. 

An early successful model of the vortex shed from the leading edge of a delta wing is 
described in Reference 7. In this model, the feeding sheet from the wing leading edge is a 
straight vortex sheet ending in a concentrated vortex. The effect of the feeding sheet on the 
crossflow plane is neglected; however, the results are very good for low aspect ratio wings. 
A more detailed treatment of slender delta wings is presented in Reference 8. In this 
method, the effect of the feeding sheet is considered, but a conical flow assumption is 
required which limits the application to low aspect ratios. A more complex vortex model 
consisting of a cloud of discrete vortices is described in Reference 9 . This model is more 
flexible in the shape and influence of the feeding sheet and the rolled-up vortex, but it is 
applicable only to the calculation of the force and center of pressure on thin wings and 
bodies. This approach has a number of features which make it desirable for bodies with 
chines, and it is very compatible with the smooth body vortex shedding technique. 

The next improvement to the vortex cloud approach is the representation of the wing 
with a panel method as described in Reference 6. This method has the capability to predict 
wing pressure distributions under the influence of the shed vorticity, and the capability to 
consider more complex configurations is available. However, a conical flow assumption is 
part of the method, and it is not clear if the method is applicable to very small wings or 
chines. Finally, the use of three-dimensional panel methods to represent both the wing and 
shed vortex have proved to be very accurate and applicable to a wide variety of 
configurations (Refs. 10 and 11). These methods require an iterative calculation procedure 
to converge on the proper vortex strength and position which causes the run times to be 
excessive in many cases. 

The next level of complexity beyond the methods discussed above involves solutions of 
the Euler and Navier-Stokes equations. There is no doubt that the state of the art in 
computational fluid dynamics is changing dramatically at this time, and the capability to 
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handle the forebody with chme is near or even available for research purposes. However, 
for preliminary design analysis, the CFD methods are still much too expensive and difficult 
to use. These methods were not considered for the current effort. 

Before a prediction method is developed to calculate the vortex shedding from 
forebodies with chines, a fundamental decision must be made as to whether the 
configuration is a noncircular body with sharp edges or a body with a low aspect ratio 
highly swept wing. As discussed later, this distinction affects how the separation vortex 
wake is perceived and developed. In this effort, in light of the previous work on noncircular 
bodies, the chine configuration will be considered a noncircular body with a sharp edge; 
however, for certain purposes (e.g., the starting problem) it will be permitted to have certain 
characteristics of a wing or leading-edge extension. 

General Approach 

The development of an engineering method to predict the pressure distributions on 
arbitrary missile or aircraft forebodies in subsonic flow at high angles of incidence requires 
the use of a number of individual prediction techniques. In the remainder of this section, 
the individual methods are described briefly, and the section concludes with a description 
of the complete calculation procedure. 

The following analyses are extensions of the work described in Reference 1; therefore, 
much of it is repeated from that reference and included here for completeness. The specific 
modifications required by the chines will be noted as will the minor changes in the method 
which will affect the prediction of separation from smooth bodies. 

Conformal Mapping 

The crossflow plane approach applied to arbitrary missile and aircraft forebodies 
results in a noncircular cross section shape in the presence of a uniform crossflow velocity 
and free vortices in each plane normal to the body axis. The procedure used to handle the 
noncircular shapes is to determine a conformal transformation for mapping every point on 
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or outside the arbitrary body to a corresponding point on or outside a circular body. The 
two-dimensional potential flow solution around a circular shape in the presence of a 
uniform flow and external vortices is well known and has been documented numerous 
places in the literature (Refs. 12 and 13). Thus, the procedure is to obtain the potential 
solution for the circular body and transform it to the noncircular body. Conformal 
transformations used are of two distinct types, analytical and numerical. 

Analytic transformation .- For very simple shapes like an ellipse [Fig. 2(a)], the 
transformation to the circle can be carried out analytically as described in Reference 13. For 
example, 


a “ v 



( 1 ) 


where 


a = y + iz 


( 2 ) 


in the real plane, and 


v = r + iX 


(3) 


in the circle plane. Two derivatives of the transformation which are required for velocity 
transformations discussed in a later section follow 


da 

dv 


1 




(4a) 


dv 

dx 


^ db da 

b -t— - a -r— 
dx dx 



2 2 

2 „ „ 


a - b 

2 ,2 


v + 

- a + b 


4 v 



1/2 


(4b) 
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A simple wing-body configuration can be considered in the same manner using 
transformations also available in Reference 13. The transformation for an ellipse is included 
in the code; therefore, it is shewn above for illustrative purposes. 

Numerical transformatio n.- For complex noncircular shapes, the transformation cannot 
be carried out analytically and a numerical transformation is required. An appropriate 
numerical transformation procedure which maps the region outside any polygon to the 
outside of a circle is described in detail in References 14 and 15. A brief summary of the 
conformal mapping procedure follows. 

The sequence of events in the numerical mapping is shown in Figure 2(b). The 
arbitrary cross section shape of the body in the a-plane is required to have a vertical plane 
of symmetry. The transformation of interest will map the region on and outside the body in 
the a-plane to the region on and outside a circle in the v-plane. The first step is a rotation to 
the 5 -plane so that the cross section is symmetric about the real axis. Thus, 

$ = i a = <; + i n (5) 

A mapping that transforms the outside of the body in the 5 -plane to the outside of the unit 
circle is 


5 


m 

^ n (v 
V r-1 


v ) 
r 


7 /* 


dv 


(6) 


where y r are the exterior angles of the m segments of the body cross section. In 
Equation (6), II denotes a product series. For a closed body 


m 

E 7 “ 2 7r ( 2 ) 

r-1 r 

The transformation implimented is an approximation of Equation (6), and it has the form 
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a 


( 8 ) 


= -l 


K 

+ I 
k-0 


k+1 * 

r A i 

o k 


where the A k coefficients are obtained through an iterative scheme described in Reference 
14, and r o is the radius of the circle in the v-plane. The derivatives of the transformation 
required for velocity calculations described in a later section are 


da 

dv 


-l 


1 


K 
- L 
k-0 



k+1 

v 


(9a) 


dv 

dx 


K (k + 1) r A ,1 
' o k 

dr 

o 

4. 

K r cLA, 1 

o k 

2 -. v 

k-0 v J 

dx 

k dx 

k-0 v J 


K r k+1 kA, 

1 - z — k 

k-0 


k+1 


(9b) 


The above numerical mapping procedure has been applied to a wide range of general 
cross section shapes with good success. The mapping works best for shapes which can be 
considered polygonal, but it does not always handle bodies with negative exterior angles 
well. More details regarding its use are presented in the description of the code in a later 

section. 


Body Model 

A three-dimensional representation of the missile volume is needed for purposes of 
predicting the absolute pressure coefficient on the surface (Ref. 2). Since the model must be 
a computationally efficient means of representing both circular and noncircular bodies in 
compressible flow up to the critical speed, a method using discrete singularities on the body 
axis is described in this section. As described in Reference 16, a panel method can be used to 
represent the body surface; however, increased computational requirements make the 
present approach more desirable for an engineering prediction method in which numerous 
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calculations may be required for preliminary design and analysis. Though the panel 
method has advantages for noncircular cross sections, an approximate method using 
axisymmetric singularities with the conformal transformations described previously has 
proved to be a reasonable approach for the present method in both surface pressure 
calculations and vortex tracking. 

Circular bodies .- The volume of an axisymmetric body is well represented by a series of 
point sources and sinks distributed on the axis. A number of models in varying degrees of 
complexity are available for this task; for example, see References 17-21. For use in 
VTXCHN, the same discrete source/sink model (Refs. 1 and 20) used for the code VTXCLD 
was selected for its accuracy, economy, and reliabilty. A new method described in 
Reference 21 appears to have some promise in modeling axisymmetric bodies, and it should 
be considered as a possible improvement to the current approach. 

Given a series of K point sources and sinks distributed on the missile axis, the induced 
axial and radial velocities at a point (x,r) are 


K 
u 

r 


V cosa 
oo C 

k-1 

K 


V CO S a 
oo C 

k=l 






where 


Q k = 


4 7tL V cosa 


( 10 ) 


(id 


( 12 ) 


is the dimensionless source strength. Q,; and x k are the source strength and axial locations, 
respectively, of the k-th source. The body surface slope at the j-th body point obtained from 
Equations (10) and (11) is 
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dr 

eg 

dx 


j 


j 


V COSa 

CD C 


1 + 


u 

r . 

2 

V cost* 

oo C 


for j-1. 


(K-3) 


(13) 


Equation (13) produces a set of K-3 linear equations in the unknown source strengths. For a 
closed body, imposing the condition of the sum of the source strengths to be zero gives the 
relation 


K , 

L Q k " 0 (14) 

k-1 


The remaining two conditions are the imposition of stagnation points at the nose and tail of 
the body. These conditions from Equation (10) are 


K 

£ 

/ 

Q k 



r x 

2 

k=l 

k 



L 



K 

L 

k-1 



(16) 


Given the source positions on the missile axis, Equations (13) through (16) comprise the total 
set of equations to solve for the source strengths. The predicted body shape from the stream 

function is 



( 17 ) 
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which must be solved by deration. This procedure is used in the code to calculate the body 
shape for comparison with the input shape. 

The above method has proved successful in modeling a variety of axisymmetric missile 
bodies; however, some care is required because of the ill-conditioned matrix. Best results 
are obtained if the discrete 1 sources are spaced at intervals of 60% of the local radius and if 
the surface slope description of the body is smooth without discontinuities. Boundary 
conditions are satisfied at points midway between the source locations. The prediction of an 
appropriate source distribution is an automated part of VTXCLD, and the user must only 
specify body geometry. 

Noncircular bodies .- An appropriate body model for missiles with noncircular cross 
sections is a surface panel method similar to that described in References 3 and 16; however, 
the use of a panel model adds significantly to the cost of each computation. For this reason 
an alternate approach was selected for use with noncircular bodies. 

The noncircular body is replaced with an equivalent axisymmetric body having the 
same cross sectional area distribution as the actual body. There are approximations 
involved with this model since the induced u-velocity due to noncircular effects is obtained 
from two-dimensional considerations as described in a later section. This approach is based 
on a high angle of attack version of slender body theory described in Reference 22. For 
example, elliptic cross section bodies at a = 0° were modeled with both a panel method and 
the axis singularity method with correction for noncircular effects, and the surface pressure 
results are nearly identical. These results are presented in a later section. 

Compressibility effec ts.- The selection of a compressibility correction scheme must take 
into consideration the configurations of interest and the calculation procedure. For 
example, vortex shedding and tracking is very dependent on cross sectional shape, and 
numerical transformations add significantly to computation time; therefore, it is 
advantageous that the compressibility transformation used have no effect on the cross 
sectional shape. With this guideline and based on similar experiences and requirements 
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(Ref. 20), a transformation which modifies only the axial coordinate is needed. The Gothert 
Rule, described in Volume I of Reference 23 and others, is the choice for the compressibility 
transformation. A brief description of the method included in VTXCHN follows. 

The transformation from the compressible (x,y,z) coordinate system to the 
incompressible (x',y',z') system is 


x' = x/0 

y' = y 

z' = z 


where 


P - I 1 - M 2 
00 

As a result of Equation (18), the body slope and angle of attack become 


— - p — 
dx' p dx 


a = pa 

and the velocity fields are related by 

U = u'/p 2 

V = v7 P 
w = w' / P 


(18) 


(19) 


( 20 ) 


( 21 ) 


(22) 
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The vortex shedding scheme in VTXCHN requires that the compressibility correction 
be applied in a manner slightly different from the usual. At the initiation of the calculation, 
the actual missile body is transformed to the incompressible shape by the above stretching 
procedure, and the source distribution is obtained as described previously. The modified 
incompressible flow conditions are used to calculate the velocity field associated with the 
body. This velocity field is transformed back to the compressible domain (Eq. (22)) so that it 
is available for use in vortex tracking and compressible pressure coefficient calculations. 
The compressible pressure coefficient is used to locate separation and define the shed 
vorticity. Since the compressibility transformation has no direct effect on the cross sectional 
shape or the discrete wake vortices, the separation and vortex tracking calculations take 
place in the compressible domain. No compressibility effects are considered for the vortex- 
induced velocities. 


Vortex Shedding Model 

The body vortex shedding model for forebodies with chines is significantly different 
from the model for smooth bodies described in References 1 and 2. The major difference is 
in the formation of the discrete vortices themselves. Since the chine is the origin of 
separation, it is no longer necessary to predict the separation location; however, there are 
other problems associated with the specification of the vortex characteristics that are 
different for chines. These will be described in this section. 

In the development of the prediction method for forebodies with chines, some 
improvements to the vortex shedding method for smooth bodies were determined. 
VTXCHN can be used for smooth bodies, but these improvements are noted, and the 
recommendation is made that they be included in the code VTXCLD (Ref. 1). 

Equations of motion .- The equations of motion of a shed nose vortex in the presence of 
other free vortices in the vicinity of a body in a uniform stream follow. In the circle (v) 
plane, the position of a vortex, T m , is 

v - r + iX (23) 
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and an image of f m is located at 


2 

r 



m 


to satisfy the flow tangency condition on the body surface. In the real plane, the position of 
the vortex r is 

m 



y + iz 
m m 


The complex potential in the real plane is 


(25) 


W(cr) = if 


(26) 


and the corresponding velocity at T m is 


v - iw 
m m 


dW (a) 
m 

da 


d_ 

dv 


W (a) 
m 


dv 
, da 


a* a 


V“V 


m 

m 


(27) 


The complex potential of r m is not included in Equation (27) to avoid the singularity at that 
point. The derivative of the transformation required in the above equations is obtained 
from Equations (4) or (9). 


The differential equations of motion for T m in complex form are 


da v - iw 

m mm 

dx V cosa + u 

oo c 

where 


(28) 


a 


m 


- iz„ 


(29) 
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Therefore, the two equations which must be integrated along the body length to determine 
the trajectory of f m are 



dx V cosa + u 
co c 


(30a) 


and 



dx V cosa + u 
oo c 


(30b) 


where u is the axial perturbation velocity in the flow field. Details of the calculation of the 
required velocity components are discussed in the following section. 


There are a pair of equations like (30) for each vortex in the field. As new vortices are 
shed, the total number of equations to solve increases by two for each added vortex. These 
differential equations are solved numerically using a method which automatically adjusts 
the step size to provide the specified accuracy. 

Velocity Field .- The velocity components at all points in the flow field are needed for 
pressure calculation and vortex tracking. The following summarizes the procedure for 
determining the u, v, w-componen ts of velocity from all singularities in the flow field. For 
purposes of this section, the velocity components required for vortex tracking are shown to 
illustrate the procedure. 


The complex potential in the crossflow (circle) plane consists of the following 
components. 


Crossflow due to a: 


W^a) = -ivV^sina 


(31) 
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Crossflow due to 0: 


W 2 (a) = -vV^sin/3 (32) 

Cylinder in a flow (two-dimensional doublet): 

w 3 ( a) = i (r^ / vJV^sinor (33) 

Cylinder in p flow (two-dimensional doublet): 

w 4 (o) = - (r^ /v)V m sinp (34) 

Discrete vortices outside a circle: 

N r 

W 5 (a) “ _i Z 21 ln (v ' v n > (35) 

n-1 

Images of discrete vortices: 

NT 2 

W 6 (a) “ 1 Z 27 ln ( v “ ~ I (36) 

n=l 1 v ’ 

n 

Note that the image at the center of the circle as required by the circle theorem is 

omitted from the potential (Ref. 24). When a vortex is shed from a cylinder, it must leave an 

equal and opposite circulation on the cylinder. This is contrary to the situation when an 
external vortex is brought to the cylinder from infinity, in which case, the center image 
vortex is required. When a vortex or cloud of vortices breaks from the feeding sheet in an 
asymmetric condition, then it may appear to downstream planes that these vortices did not 
originate on the body, and the center vortex is required as part of the solution. In most 
practical missile situations when the shedding is symmetric, the above point is academic to 
the method. 
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The next term is that due to a two-dimensional source representing a growing cylinder 
in the crossflow plane. It is required for the incremental noncircular perturbation effects on 
the u-velocity to be described later. Note that the axisyrnmetric body contribution to all 
velocity components is included through the three-dimensional source/ sink distribution 
described previously in Equa tions (10) and (11). 

Expanding cylinder (two-dimensional source): 

dr 

W (< 7 ) =* r lnv V cosa (37) 

7 eq dx 

The velocity component in the crossflow plane are obtained from the derivative of the 
complex potential as shown in Equation (27), The total velocity at T m in the crossflow 
plane is written as 


v - iw 

— = G +G. + G + G + G + G (38) 

V a B n m T r 

oo 

where each term in Equation (38) represents a specific velocity component in the d-plane 
corresponding to the appropriate complex potential. The first term, from Equations (31) and 
(33), represents the uniform flow due to angle of attack. 
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The second term, from Equations (32) and (34), represents the uniform flow due to 
angle of yaw. 
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In compressible flow, the velocity components from Equations (39) and (40) used in the 
surface pressure calculations include effects of the compressibility transformation described 
in a previous section. 

The next term, from Equations (35) and (36), represents the influence of all vortices and 
their images, with the exception of T m . When velocity components are required elsewhere 
in the field, T m is included in the calculation. 
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The next term, from Equation (36), is due to the image of r . This term is required 
only for vortex tracking purposes. 
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The fifth term in Equation (38) represents the potential of T m in the a~plane (Routh's 
Theorem, Ref. 24) and is written as 
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(43) 


The last term in Equation (38) corresponds to the velocity components induced by the 
three-dimensional source singularities representing the volume effects from Equation (11). 
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oo 


( 44 ) 


where v r and w r are appropriate components of v R . 



It has been demonstraled with panel methods that the u-velocity perturbation due to 
noncircular effects is important in the calculation of surface pressures and flow field 
velocities (Ref. 16). Since all the singularities in the flow model are two-dimensional, with 
the exception of the axisymmetric body volume model, it is necessary to calculate the u- 
velocity components from the complex potentials defined in Equations (31) through (34) and 
(37). The u-velocity contribution of the vortex field is neglected in the flow field, but it is 
included in the surface pressure calculation as described tn a later section. The following is 
a description of the calculation of the various components of the u-velocity at any point in 
the flow field. 


The axial velocity, given the complex potential, is 
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From Equations (31) through (34), the noncircular contributions to Equation (45) are 
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Since the body source singularities given by Equation (12) are three-dimensional, they 
contribute an induced axial velocity, u r , given by Equation (10). However, as noted 
previously, these axial velocities are axisymmetric, and they exhibit no direct effect of the 
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noncircular shape. This deficiency can be corrected using the techniques of high angle-of- 
attack slender-body theory described in Reference 22. An approximation to the 
perturbation u-velocity due to the growing noncircular shape is described below. 

The equivalence rule described in Section 6-4 of Reference 25 states that (1) the flow far 
away from a general slender body becomes axisymmetric and equal to the flow around an 
equivalent axisymmetric body, and (2) near the slender general body, the flow is different 
from that around the equivalent axisymmetric body by a two-dimensional component that 
is required to satisfy the tangency boundary condition. Or, the difference between the 
velocity potential for the noncircular body and that for the equivalent body of revolution is 
equal to the increment of the velocity potential due to the noncircular effect. Thus, 
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(50) 


Applying Equation (45) to Equation (37) for both terms on the left side of Equation (50), 
remembering that dv/dx = 0 for circular cross sections, the contribution due to the growing 
noncircular forebody becomes 


dW ? 

dx 


= V cosa 

r dr 

eq eq 

dv 

oo 

An 

v dx 

dx 


(51) 


This velocity represents the increment in the axial velocity caused by the noncircular shape 
of the body. 


Finally, the u-velocity in the flow field of a noncircular body is calculated by including 
Equations (46) through (49) and (51) in Equation (45). This approach, though slightly 
different from that presented in Reference 22, produces results that are in excellent 
agreement with high angle-of-attack slender-body theory as applied to ellipsoids. 

Surface pressure distribution .- The surface pressure distribution on the body is 
required to calculate the forces on the body and the separation points. The surface pressure 
coefficient is determined from the Bernoulli equation in the form 
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where U is the total velocity (including VJ, velocity components from Equation (38), and u 
velocity components from Equations (46) through (49) and (51) at a point on the body. 


The last term in Equation (54) represents the axial velocity missing from the two- 
dimensional vortices representing the shed vortex field. This "unsteady term in 
Equation (54), evaluated on the body surface, is 
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The contribution of the discrete vortices in the flow field, using Equation (35) and (36), 
becomes 
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where 


There are several options available for the use of the information defined in Equation 
(56) and in the previous section. Practical aspects of the calculation of pressure coefficients, 
velocity fields, and vortex trajectories dictate the manner in which these velocities are most 
efficiently considered. This is discussed in detail in the section on Calculation Procedures to 

follow. 

Separated wake .- The separated wake on the lee side of the body is represented by a 
large number of discrete vortices, each vortex originating from separation locations at each 
axial marching step in the calculation. The major portion of the lee side vortex wake has its 
origin at the primary separation points on each side of the body. The remainder of the wake 
originates from the secondary separation points located in the reverse flow region on the lee 
side of the body. Both of these points are illustrated in the sketch of a typical crossflow 
plane of an elliptic cross section body shown in Figure 3. The mechanics of the calculation 
of the individual vortices for both smooth forebodies and forebodies with chines are 
described to illustrate the differences in the prediction methods. 
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For smooth bodies. References 1 and 2, the separation causing the vortex shedding into 
the lee-side wake is of a type traditionally associated with boundary layer separation. 
Prediction of the separation locations on smooth bodies has been presented in the above 
references, but a brief descrip tion is included here for completeness. 

The predicted pressure distribution for the primary flow in the crossflow plane is 
referenced to the conditions at the minimum pressure point, and a virtual origin for the 
beginning of the boundary layer is calculated. The adverse pressure distribution 
downstream of the minimum pressure point is considered with either Stratford's laminar 
(Ref. 26) or turbulent (Ref. 27) separation criterion to determine whether or not separation 
has occurred. These criteria, based on two-dimensional, incompressible, flat plate data, are 
adjusted empirically for three-dimensional crossflow effects in Reference 2. The laminar 
separation criterion states that the laminar boundary layer separates when 
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In a turbulent boundary layer, separation occurs when 
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(Re x 10 * 0.350 sina (59) 

£ c 

The sina c in Equations (58) and (59) is the three-dimensional modification. 



If the criteria indicate a separation point, the vorticity flux across the boundary layer at 
separation is shed into a single point vortex whose strength is 


r 

v 


2 

u Ax 
e 


2 

2V cosa 

oo C 


6 


( 60 ) 


- 27 - 


assuming no slip at the wall. 6 is the vortex reduction factor described previously, and it is 
generally set equal to 0.6 for ogive-cylinder configurations and 1.0 for bodies with closed 
boattails (Ref. 1). For bodies with chines or other sharp edges initiating separation, 6 = 1. 

For smooth bodies, the initial position of the shed vortex is determined such that the 
surface velocity in the crossflow plane at the separation point is exactly canceled by the shed 
vortex and its image. When this criterion results in a vortex initial position that is too near 
to the body surface, certain numerical problems can cause difficulty in calculating the 
trajectory of this vortex. If the initial position of the vortex is nearer than five percent of the 
body radius from the surface, the vortex is generally placed five percent of the equivalent 
body radius from the surface. 

For bodies with chines or other sharp edges, the primary separation location is fixed at 
the sharp edge, thus negating the need to predict separation as described above for smooth 
bodies. The appropriate boundary condition is to assume the flow leaves the sharp edge 
smoothly; that is, the Kutta condition is satisfied at the edge. This boundary condition 
transforms to a stagnation point in the circle plane at the separation point. 

The common difficulty for all vortex cloud methods applied to sharp edges (e.g,, 
Refs. 5, 6 and 9), including the present method, lies in the solution for the shed vortex at 
each time step. Only one equation is available from the Kutta condition to solve for two 
unknowns, the vortex strength and its position. In the smooth body analysis, an external 
boundary condition provides the missing equation. This condition is the determination of 
the strength of the vortex from the vorticity in the boundary layer prior to separation as 
given by Equation (60). In some methods noted above, the analysis was developed with the 
assumption of conical flow to provide the necessary information. Several different 
approaches were examined in the present work, and these are discussed below. 

The initial approach is based on the method developed by Sacks (Ref. 9). At the first 
shedding station on the forebody and chine, there is no vorticity in the field; therefore, the 
chine edge is a singularity and it is not possible to calculate a realistic v-velocity near the 
edge. In Reference 9, the solution is started by assuming the velocity outboard at the chine 
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edge is approximately represented by v = 0.5sina, and the shed vortex is convected 
outboard a distance vAt. In this work, the pressure distribution is calculated at the initial 
station with no vorticity in the field; therefore, there is a singularity at the chine edge. The 
average v-velocity over the outer 10-percent of the chine, excluding the singularity, is used 
to locate the initial shed vortex. Placing the vortex at this position satisfies the requirement 
for one of the two equations, and the remaining equation can be used to calculate the 
strength of the vortex such that the Kutta condition at the chine edge is satisfied. With a 
vortex in the field which satisfies the Kutta condition, it can be convected with the local flow 
field during the next time step. Now, though the Kutta condition is not satisfied after the 
vortex moves, the v-velocity can be calculated at several points on the forebody near the 
chine edge. The velocities on the outer ten percent of the chine semispan are averaged to 
find a convecting velocity for the next shed vortex. As before, when the position is known it 
is possible to determine the strength by satisfying the Kutta condition at the chine edge. 

The method of convecting the shed vortex to its initial position has been the subject of 
much discussion between investigators, and there are nearly as many methods as there are 
investigators. When a sheet of vorticity is shed from the trailing edge of an airfoil in 
unsteady motion, a case is made in Reference 28 that the sheet should leave the edge tangent 
to the surface. In Reference 29, a procedure for locating the shed vortex next to the edge of a 
cambered plate is described. In this case, the vortex is located at a point in the field which 
produces the proper flow field in the vicinity of the edge. Unfortunately, this approach 
requires a significant iterative computational effort for each different geometry; therefore, it 
is not practical for an engineering method which must be available to use for changing 
geometries without large set-up times. In the current version of VTXCHN, the former 
method is used in which the shed vortex is convected to its initial position along a line 
tangent to the windward surface of the chine. 

An alternate approach in which the chine is approximated as a lifting surface is also 
available to start the solution. In this method, the chine is represented as a single panel 
lifting surface on the transformed circular body as shown in Figure 5. The single panel is 
oriented along the mean line of the chine to approximate the local incidence to the free- 
stream velocity. The no-: low boundary condition is satisfied at a single point on the panel, 
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and the strength of the vortex at the panel edge is the result. With the strength known, 
there is a single point along the projection of the panel used to model the chine at which the 
vortex can be located to satisfy the Kutta condition at the chine edge. This provides the 

necessary starting solution so that the remainder of the calculation can be carried out as 
described above. 

The mechanics of the panel approximation are very simple. After some investigation, 
the best location to satisfy the boundary condition is at half the chine semispan. It is also 
best to use this approximation only for the starting solution, even though it can be applied 
at any number of axial stations. This latter limitation is imposed because the strength of the 
vortex determined from the panel solution is significantly weaker than the vortex shed by 
the previously described method. For example, in Figure 6, the build-up of shed vortex 
strength from the two methods is compared for a typical chine forebody configuration at 
a = 20* . Notice that the strength differential caused by the initial shed vortex specified by 
the panel approximation is never recovered, but the difference over the length of the 
forebody is minimal and has almost no effect on the forebody loads. 

The chine shedding solution was further studied by specifying the strength of the shed 
vortex in a manner similar to that used for smooth bodies. The velocity distribution tangent 
to the lower surface of the chine near the outer ten-percent of the semispan was calculated 
and averaged to approximate the velocity near the chine edge. Assuming that separation 
occurs at the chine edge, but also assuming that the vorticity in the boundary layer at the 
chine edge determines the strength of the shed vortex as on a smooth body, a vortex can be 
placed in the outer flow to satisfy the Kutta condition at the edge as described above. A 
preliminary calculation with this method resulted in a total vortex strength very similar to 
the other methods described, but the vortex shedding rate seemed to be erratic. This may be 
caused by the difficulty in calculating a smooth velocity distribution necessary to determine 
vortex strength in a region where the local velocities can change quickly as the discrete 
vortices move in and out of influence of the chine edge. There was also difficulty satisfying 
the Kutta condition at the chine edge; therefore, this method was abandoned. 
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The calculation of secondary separation on the lee side of a forebody with a chine is 
carried out in the same manner as described above for primary separation on smooth 
bodies. It is necessary that a reverse flow region exist on the lee side of the body and that a 
second minimum pressure point be found in this region. For purposes of this analysis, the 
reverse flow is assumed to be laminar from the lee side stagnation point to the secondary 
separation point, and Stratford's laminar criterion is used to locate secondary separation. 
Laminar separation in the reverse flow region is expected because of the low velocities on 
this portion of the body. The vortex released into the flow at the secondary separation point 
has the opposite sign of the primary vortex and is generally weaker in strength. 

In theory, secondary separation is a straightforward concept; however, in practice, it is 
very difficult to implement. This has been noted by other investigators (Ref. 29) for various 
body shapes. The computational difficulties are caused by both the separation and the 
tracking portions of the method. The separation problems are caused by the variations in 
the pressure distributions on the lee side of the chine caused by the influence of the primary 
chine vortex. It is not unusual that a second minimum pressure is not observed at some 
axial stations; as a consequence, a secondary separation point is not predicted and a "hole" is 
left in the secondary vortex sheet. This problem can be overcome by forcing secondary 
separation, a common procedure in vortex cloud methods. The tracking problems are 
caused by the large primary vortex capturing the smaller secondary vortices and pulling 
them away from the chine surface. When this happens, the induced effects of the secondary 
separation disappear, and there is no advantage to including these additional vortices in the 
solution. In addition, though not a primary consideration in this discussion, the secondary 
vortices have a significant influence on the computation time since each vortex adds two 
more differential equations to the set to be solved. 

Forces and momen ts.- The forces and moments on the body are computed by 
integration of the pressure* distribution around the body. At a specified station on the body, 
the normal-force coefficient on a Ax length of the body is 
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where T is the distance from the axis of the body to the body surface and 0' is the local slope 
of the body in the crossflow plane. This is illustrated in the sketch in Figure 4. For circular 
bodies, r' = r Q and 0' = 0. The total normal force coefficient on the body is 
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and the pitching-moment coefficient is 
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The center of pressure of the normal force is 
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Similarly, the side-force coefficient on a Ax length of the body is 
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The total side-force coefficient on the body is 
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and the yawing-moment coefficient is 
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The center of pressure of the side force is 
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A noncircular body at an arbitrary roll angle may experience a rolling moment caused 
by the nonsymmetry of the loading around the body. The total rolling moment on the body 
is calculated by summing the moments of the individual components of normal force and 
side force around each cross section and integrating over the body length. The rolling 
moment coefficient is calculated as 
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Though the primary purpose of the vortex shedding method is not to predict drag or 
axial force on the missile, a procedure to estimate both pressure drag (excluding base drag) 
and skin friction is included. The pressure contribution to the axial-force coefficient is 
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The local skin-friction coefficient, based on the assumption of a 1 /7th power law velocity 
profile in the boundary layer, is 
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which produces a drag coefficient due to friction 
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At high angles of incidence, streamlines around the body are inclined at approximately 


a s = tan _1 (2tana ) 


(73) 


therefore, the axial component of the friction drag is 
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The total axial-force coefficient is the sum of Equations (70) and (74). 


Vortex core. - The diffusion core model (Ref. 2) for the point vortex- induced velocities 
removes the singularity at the vortex origin and effectively reduces the velocities near the 
vortex. The tangential velocity induced by a single point vortex is written as 
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where r is the distance from the vortex to the field point and x is a measure of the age of the 
vortex. The induced velocity from this core model is illustrated in the following sketch. 


- 34 - 



Sketch - Vortex-Induced Velocity 

The vortex core model represented by Equation (75) has received considerable 
attention in the context of body vortex-induced effects, and it has a number of 
shortcomings. Since the exponential term is a function of r, the flow medium (v) and the 
age of the vortex (x), the core is constantly changing size as the vortex moves through the 
field. Under certain conditions, the core radius, denoted as r c in the sketch, can become 
very small and the induced velocity becomes unrealistically large. In an attempt to further 
modify the core model to keep the induced velocities within physically realistic limits, the 
following modification was made in Reference 3. The location of the maximum induced 
velocity, r c in the sketch, is fixed at a specific radius to be selected by the user. Given a core 
radius, the vortex-induced velocity is 
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The core model described by Equation (76) is included in the version of the code 
described in this report, and in a later section, guidelines for the selection of an appropriate 
core radius are presented. Results indicate that this simple core model provides adequate 
smoothing for the necessary velocity calculations. 


Other investigations of discrete vortex models must use a core model of some type. 
Although there are many different core models, ail serve the same purpose as those 
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described above, and nearly all are directed at eliminating the singularity and reducing the 
maximum induced velocity The code VTXCHN described in this report is easily modified 
to incorporate another core model if the user desires. 


Calculation Procedures 

As a part of the chine investigation described herein and based upon the experience of 
the authors in working with vortex cloud methods, several interesting changes in the 
calculation procedures have been studied. These changes were identified during a study of 
computational efficiency and accuracy in the anticipation that both could be increased. The 
following paragraphs discuss the results of this study. 

Unsteady pressure term .- Practical aspects of calculating the pressure distribution on 
chined forebodies for preliminary design and analysis dictates the requirement for a fast 
computational method. As has been described by a number of investigators, discrete vortex 
methods can be time consuming and computationally inefficient, particularly when the 
number of vortices is large. For example, prediction of the pressure distribution from 
Equation (56) can require a significant amount of computer time on bodies with a large 
number of vortices simply because each vortex contributes to this term at each point on the 
body. When a noncircular body is being considered, the computation time is increased 
further because each vortex must be mapped numerically to the circle plane. Computation 
time notwithstanding, it has also been noted by the authors that this unsteady pressure term 
can cause a number of other problems. For example, when an individual discrete vortex is 
very near the surface of the body, it can have a very large local effect on the pressure 
distribution both through the induced velocity effect and the unsteady effect. The unsteady 
pressure term can be large even when an individual vortex in the cloud has occasion to 
move rapidly in the field because of interaction with another vortex. In both these cases, the 
local pressure coefficient may have a saw-tooth distribution and can exhibit extreme jumps 
or spikes. The goal of this study was to find a means to smooth the predicted pressure 
distribution and achieve faster computational times with no significant change in accuracy. 
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Since the major contributor to the computation time of the unsteady pressure term. 
Equation (56), is the summation over all the vortices in the field for each point on the body, 
this is where the study began to increase computational efficiency. Two approaches were 
examined. In the first, all shed vorticity was concentrated into a single vortex at the 
centroid of the cloud for purposes of evaluating the vortex contribution of Equation (56). In 
the second approach, the feeding sheet was represented by a number of discrete vortices 
and the remainder of the cloud was concentrated into a single vortex located at the centroid 
of the cloud. This latter approach is illustrated in Figure 7. 

Before any evaluation of computational savings of the modifications, it is important to 
understand the effect on the accuracy of predicted results. To this end, comparisons of 
measured and predicted pressure distributions on an axisymmetric ogive-cylinder body 
(Ref. 30) were used to evaluate the calculation procedures. An axisymmetric body without 
chines was selected to avoid any special problems created by the noncircular shape and the 
associated mapping and to take advantage of well-tested experimental data in which 
confidence exists. 

A series of comparisons at a - 20° is shown in Figure 8 where the circumferential 
pressure distribution at a number of axial stations along the body are considered. For 
purposes of this evaluation, three different cases are shown. The first, denoted by NVPHI=- 
1, corresponds to the traditional calculation procedure in which all vortices in the cloud are 
included in the summation in Equation (56). The second, denoted by NVPHI=0, represents 
the condition in which all shed vortices at each axial station are collected into a single vortex 
at the centroid of the vortex field. Finally, the last case, denoted by NVPHI=5, represents a 
combination of the previous two cases in which the feeding sheet consists of five vortices 
and the remainder of the cloud is collapsed into a single vortex. Remember that the 
individual vortices in the cloud are collapsed into their centroid only for the evaluation of 
Equation (56), and for all other computations, the entire cloud of discrete vortices is used. 

In Figure 8(a), there is no difference in the predicted pressure distribution near the 
nose. In Figure 8(b), there is a difference in the pressure distribution between the traditional 
results (NVPHI=-1) and those results obtained when all the discrete vortices are collapsed 
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into a single vortex at the centroid (NVPHI=0). There is no effect of keeping the feeding 
sheet and cloud separate because the total number of vortices in each field is less than seven; 
therefore, this method and the traditional method are identical at this point. 

As the number of vortices shed from the body increases, the difference between the 
traditional results and the alternative methods becomes larger, and as seen in Figures 8 (c, d, 
and e), collapsing all the vortices into a single vortex produces significantly inferior pressure 
results. Keeping the feeding sheet intact and collapsing the cloud vortices is in better 
agreement with the original method, but it appears that, as the number of vortices becomes 
large, the agreement between the methods begins to deteriorate. It is apparent from these 
comparisons that the success of the vortex cloud method in predicting accurate pressure 
distributions is due to the distribution of the discrete vortices throughout the flow field. 

As a consequence of these comparisons, the decision was made to keep the traditional 
approach of separate and distinct discrete vortices for the entire pressure calculation. The 
feeding sheet is obviously an important feature of the flow model, and the minimal savings 
in computational efficiency observed is not worth the loss of accuracy in the predictions. 
The options to use these modifications have been retained in the code in the event the user 
wishes to test the modified method for other configurations. 

Vortex tracking .- Calculation of the motion of the discrete vortices after they are shed 
from the body is a key component of the prediction method. The location of the vortices 
and the cloud formed influences the strength and subsequent positions of later vortices shed 
from the body, and ultimately, the tracked positions of the vortices determines the induced 
loading on the body. As noted by other investigators (Ref. 27), it is common for tracking 
problems to arise during the normal calculation of the motion of the individual vortices. 
The usual nature of these problems is a vortex inside the body, and the usual solution is to 
simply eliminate the offending vortex from the flow field. In the code presented herein and 
in Reference 1, the solution is to replace the vortex back into the flow field outside the body 
and allow it to continue as part of the cloud. This procedure has been reasonably successful 
on axisymmetric bodies, though the problem occurs so infrequently in most configurations 
examined by the authors that it is difficult to test the validity of this approximate solution. 
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This procedure did not prove very successful in the case of the chine configurations; 
therefore, some work was directed at understanding the mechanics of the tracking problem 
in the anticipation that it could be eliminated. 

The equations of motion of the vortices are shown in Equation (30). The motion of a 
vortex is determined by integrating these equations, along with those for the other vortices, 
from one axial station to the next, a distance dx. For most purposes, dx is a constant over 
the entire length of the body; although, as discussed in Reference 3, dx can be a variable 
length. When the tracking procedure was originally developed, the u-velocity in the 
denominator was included to correct for the effect of a growing body near the nose; that is, 
for axisymmetric bodies, the u-velocity was determined by the source singularities 
representing the body (Eq. (10)). This velocity was essential for the successful tracking of 
vortices shed near the nose, and, upon a brief analysis, it becomes obvious why this is the 
case. Equation (30) represents to first order the slope of the vortex filament between the two 
axial stations. If a vortex is very near the body surface, the slope must be the same as the 
body slope or greater, or the vortex will be tracked inside the body; therefore, near the body 
nose the u-velocity from the sources (which is opposite in sign from V^cosa^ ) is necessary 
to decrease the magnitude of the denominator and increase the slope of the filament, thus 
keeping the vortex from penetrating the body surface during tracking. 

For noncircular bodies, the effect of the body shape on the v- and w-velocities in 
Equations (30a) and (30b) is included by means of the conformal transformations described 
previously. In the original work (Refs. 1 and 2), the u-velocity was included as an 
axisymmetric effect from the equivalent body, and there was no noncircular effect on the u- 
velocity. In this extension of the work of Reference 1, the approximate contribution of the 
noncircular shape to the u- velocity is determined from Equation (45). When the full u- 
velocity, described in the section containing Equation (45), was included in the tracking, the 
vortices tracked inside the body on a regular basis. The difficulty appears to be a subtle 
violation of one of the basic assumptions inherent in the vortex cloud procedure. 

The unsteady analogy used in the vortex cloud method of modeling the lee-side 
vorticity assumes that all the vortices are moving at the same u-velocity such that they 
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traverse the dx-distance between axial stations in the same time. When the u-velocity 
consists of a free-stream component and a volume or source component, this assumption is 
nearly correct for all vortices near the body, even when the body is noncircular and there is 
a small variation in u-velocity around the body. However, when the doublet effects are 
included, the u-velocity can vary significantly around the body and with distance from the 
body. In this case, the vortices are not all traversing axially at nearly the same velocity, and 
the unsteady analogy approximation is not strictly correct. The numerical problems can be 
explained using Equation (30). On the lee side of the body, the doublet singularities 
represent the acceleration of the u-velocity on that side of the body, the denominator 
becomes very large near the body, and the slope of the vortex filament decreases until it can 
be less than the slope of the body. As a consequence, the vortex tracking calculation forces 
the vortices inside the body. 

For purposes of the prediction method in the code VTXCHN, the u-velocity used for 
tracking purposes does not include the doublet terms. These terms are included in the 
surface pressure calculation and the velocity field calculations. 


RESULTS 

For purposes of evaluating the accuracy and range of applicability of the subsonic 
vortex shedding model for forebodies with chines and the associated computer code 
VTXCHN, comparisons of measured and predicted aerodynamic characteristics are 
presented. Since the major objective is the validation of the prediction method, not all 
predicted results are in good agreement with experiment to illustrate to the user where 
problems may occur during general use of the code. Typical results from the prediction 
method follow. 


Circular Bodies 

Though not specifically part of the stated objective of the reported work, the prediction 
method is applicable to axisymmetric bodies and noncircular bodies without chines. The 
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code, VTXCHN, can be used to predict vortex shedding from forebodies alone in exactly the 
same manner as the predec essor code, VTXCLD, in Reference 1. For circular bodies, 
VTXCHN should produce results similar to VTXCLD, but because of minor improvements 
in the separation calculation, the results may not be identical. For practical engineering 
purposes, the results from the two codes are so nearly the same, additional comparisons 
with experiment for circular bodies will not be included in this report. 

Measured and predicted pressure distributions on an ogive-cylinder model (Ref. 30) in 
subsonic flow at an angle of attack of 20 degrees has already been presented in Figure 8. 
The configuration has a 3-caliber ogive nose and a 7.7-cahber cylindrical afterbody, and as 
described in Reference 1, the model is represented by 61 point sources and sinks on the axis, 
and the ogive nose is used to close the model at the base. Laminar separation is used for all 
predictions for this model; however, there is some indication of both laminar and turbulent 
separation in the data. 

At a = 20 6 , a reasonable vortex field is developed on the lee side and has a significant 
effect on surface pressure distributions on most of the body. Some of the roughness in the 
predicted results is caused by individual vortices moving too near the body surface during 
the trajectory calculation. The vortex core model does tend to smooth the vortex-induced 
effects, but there can still be a large local effect. These local irregularities have a minimal 
effect on the integrated loads. 


Noncircular Bodies 

The noncircular body options (without chines) in VTXCHN have been tested for a 
number of different cross section shapes, both elliptic and arbitrary, to verify the 
transformation procedures. Since these are nearly identical to those reported in Reference 1, 
they are not repeated. It should be noted that the changes in the u-velocity caused by the 
noncircular body effects now included will cause the results from VTXCHN to be slightly 
different from those obtained from VTXCLD. Examples of the predicted differences in the 
pressure distributions on a 2:1 cross section elliptic body (Ref. 31) are shown in Figure 9. In 
Figure 9(a), the circumferential pressure distribution near the nose is shown at a = 0 to 
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eliminate the influence of all but the noncircular source contribution to the pressure. The 
dashed curve is the result from VTXCLD (Ref. 1), the solid curve is the result from a panel 
method (Ref. 16), and the dotted curve is the result from VTXCHN including the 
noncircular u-velocity modification described previously. The VTXCHN result is in very 
good agreement with the panel result, and it shows a significant improvement over the 
previous result without the u-velocity correction. 

The predicted pressures on the same configuration at a = 10° are compared in Figure 
9(b). In this case, no separation is included in any of the results to permit comparison of 
potential flow effects only. Now the angle of attack effects dominate the pressures, and the 
body volume effects are small; however, the VTXCHN results are still in good agreement 
with the panel results and show an improvement over the original VTXCLD results. 

A thorough examination of noncircular bodies without chines was beyond the scope of 
the present work; therefore, the results presented in Reference 1 were not repeated for this 
report. The user should be aware that VTXCHN will produce slightly different results than 
VTXCLD because of the modifications in the noncircular u-velocity terms. 

Forebodies With Chines 

Comparison of measured and predicted aerodynamic characteristics on typical 
forebodies with chines is the most reliable means to evaluate the strengths and weaknesses 
of the engineering method embodied in the code VTXCHN. Since pressure distributions 
provide the best assessment of the capabilities of a prediction method, the configuration and 
experimental data of Reference 32 provide a wide range of information for verification of 
the method. The forebody model with the three specific axial stations at which pressures 
are measured is shown in Figure 10. The experimental data shown in the following figures 
were measured with a 60-degree delta wing in place aft of the forebody, but it appears that 
the wing had only minimal influence at the last pressure station on the body. The effect of 
the wing is not modeled in the prediction method. 
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Measured and predicted pressure coefficients on the upper surface of the forebody are 
shown in Figure 11 for the configuration at a = 10\ Predicted pressures on the lower 
surface are shown in these figures, and the predicted vortex cloud pattern at each station is 
also shown to provide some perspective of the region of influence of the lee side vortex. In 
Figure 11(a), the measured and predicted pressure coefficients at the first forebody station, 
x * 7 inches from the nose, are compared. The predicted pressures have the correct 
characteristics on the suction side of the body, but the suction peak caused by the vortex 
seems to be spread over a larger region than that measured. The pressure on the lower 
surface is nearly constant. The predicted pressure distribution without vortex effects is 
shown as a dashed curve. At this low angle of attack, the vortex-induced effects are 
confined to the upper surface of the chine. 

In Figure 11(b), the second axial station, x * 13.5 inches from the nose, exhibits lower 
suction pressures on the upper surface, and the predicted results are in very good 
agreement with experiment except for the suction peak near the chine edge. 

The last axial station at x = 20 inches in Figure 11(c) shows good agreement between 
the measured and predicted pressure distributions on the upper surface. The vortex- 
induced effects are less at this station than at the previous station even though the vortex 
strength has increased. This is caused by the increased distance between the rolled-up 
vortex and the chine upper surface. The change in vortex shape between the two axial 
stations is determined by the angle at which the vortices leave the chine edge, and it is 
apparent that the vortices from the aft station are leaving the edge at a smaller angle. This 
causes the rolled-up vortex to be displaced outboard and upward from the chine edge, thus 
decreasing the vortex-indua?d pressure. 

Measured and predicted results on the same configuration at a = 20* are presented in 
Figure 12 for the same three axial stations. As before, the predicted pressure distribution on 
the forebody without separation is shown as the dashed curve to illustrate the induced 
effects of the vortices. As expected from the vortex fields, the influence of the chine vortex 
is concentrated on the chine and the region of the chine-body junction. 
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In general, the results at a = 20° with separation effects included are the same as 
shown in the previous figures for the first two axial stations. The agreement between 
measured and predicted pressures is very good. At the aft axial station, the predicted 
suction pressure on the lee side is significantly lower than that measured. The explanation 
is not obvious at this time, but the problem could be with the predicted strength and/or 
position of the shed vortex, or it could be due to some external influence such as the wing. 
There is not sufficient experimental data available to determine the reasons for the observed 
problems. 

When the angle of attack is increased to 30 and 40 degrees in Figures 13 and 14, 
respectively, the results become very consistent at all stations on the forebody. The 
predicted pressure distribution is in good agreement with the measurements on the fuselage 
upper surface away from the chine, but the agreement is generally poor everywhere on the 
chine. The predicted suction pressure is always less than that measured; therefore, the 
integrated normal force will be less than that measured at the higher angles of attack. The 
problem on the chine could be caused by the vortex position. It is apparent from the vortex 
fields shown on each figure that as the angle of attack increases, the vortex moves farther 
from the body thus reducing its effect. Comparison of the predicted pressures with and 
without the vortex field shows that the method represents the qualitative effects of the 
vorticity even though the actual pressure levels are not in good agreement with experiment. 

It is possible that secondary separation on the lee side of the forebody can be changing 
the vorticity distribution in the flow field sufficiently to have a measurable effect on the 
pressure distribution. Secondary separation effects were not included in these calculations. 

When the forebody at angle of attack is yawed to some nonzero sideslip flow angle, the 
separation becomes asymmetric. The assumption is made that separation still occurs at the 
chine on both sides of the forebody, but the asymmetric flow field changes the tracking and 
subsequent vortex distributions. This assumption is good for moderate yaw angles, but 
when the yaw angle approaches the angle of attack, the flow is approaching the windward 
chine at a low incidence angle, and separation occurs on the body surface rather than at the 
chine edge. 
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In Figures 15 and 16, the predicted pressure distributions on the forebody at 0 = 10' 
and a - 20° and 30”, respectively, are compared with experiment. The agreement at the 
first axial station is reasonable on the chines at both angles of attack, and it is clear that the 
vortices are not having enough influence on the local pressures on the lee side of the chine. 
Some of the difficulties with these results can be explained by studying the vortex field, also 
shown in the figures. At the first station, the vortex on the right side of the body seems to be 
spread out over a large region of the lee side of the body, and the influence of the vortex on 
the surface pressure reflects this phenomenon. The predicted vortex induces a low pressure 
over more of the upper surface than is measured, but the pressure peak near the chine edge 
for a = 30° is lower than indicated by the experiment. As for the 0 = 0’ described above, 
the pressure level at the lower angle of attack is in better agreement with the measured 
pressures. On the left side of the body the vortex is rolled up in a loose fashion and is 
concentrated away from the body surface. As a consequence, the predicted suction pressure 
peak is significantly lower than that measured. Another interesting result at this first station 
is the presence of a single vortex from the right chine on the left side of the body. 

At the second axial station shown in Figures 15(b) and 16(b), the vortex on the right 
side is much stronger than at the previous station, it is tightly rolled up, and, as a 
consequence, the suction peak is in better agreement with experiment. The vortex on the 
left side is also tightly rolled up, but it is concentrated away from the body and its influence 
on the pressure distribution is less than that measured at the higher angle. As at the 
previous station, the pressure peak on the left side is less than the measured value at both 
angles of attack. 

Finally, at the third axial station shown in Figures 15(c) and 16(c), the vortex is tightly 
rolled up near the body, and the predicted pressure peak on the chine is in good agreement 
with experiment. Though the left chine vortex is rolled up a large distance from the chine 
surface, the pressure levels are in reasonable agreement with the measurements except for a 
region on the body near the junction with the chine. 
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The discrete vortex fields shown in the above pressure figures help illustrate the extent 
of the vorticity shed from the chine, but these results do not provide much information on 
the nature of the physics of the flow in the vicinity of the forebody. To better understand 
the full effect of the shed vorticity, flow field calculations provide a direct visualization of 
the velocity field. A representative set of velocity vectors in the flow field adjacent to the 
forebody at a = 20° is shown in Figure 17. The accuracy of the flow vectors inside the 
vortex or very near individual vortices will be in question because of the influence of the 
discrete vortices, but the overall character of the flow field is correct. These vector plots are 
most useful in visualizing the region of influence of the vortex, particularly if there is 
concern about vortex-induced loads on other components of the airframe. 

If there is interest in estimating these vortex-induced loads, the strength of the shed 
vorticity is needed. This detailed information is also available from the prediction method. 
The strength of the vortex shed from the forebody at various angles of attack is shown in 
Figure 18 to illustrate this capability. 


PROGRAM VTXCHN 
General Description 

This section provides a general description of program VTXCHN and its various 
subroutines. The code is written in modular form so that, as improvements in flow models 
become available or as other modifications are required, they are easy to incorporate. The 
flow of the calculation is described in this report, and detailed comments are provided 
throughout the code to assist the user in understanding the order of calculation. 

The computer code consists of the program VTXCHN and 45 subroutines. The overall 
flow map of the program is shown in Figure 19 where the general relationship between the 
subroutines and external references can be seen. Communication between the program 
modules is handled primarily by named common blocks. A cross reference table for the 
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calling relationship between the program subroutines is shown in Figure 20. A similar table 
for the named common blocks is shown in Figure 21. 

The program is written in standard FORTRAN. Execution on a VAX 11/750 can vary 
from 10 seconds to more than 45 minutes depending on various factors such as the 
geometry, the integration interval (DX), the number of shed vortices, the use of flow 
symmetry, and the flow conditions. 

Subroutine Description 

This section briefly describes the main program V rXCFIN and its various subroutines. 

VTXCHN This is the main program for calculating the forces and moments on a 
slender body in a steady flow condition including the nonlinear effects of 
lee side separation vorticity. 

AMAP Transforms a specified point in the circle plane to its corresponding point 
in the body plane. 

ASUM Calculates the velocity term due to transformation for noncircular bodies. 

BMAP Sets up a table of points on the actual body corresponding to specified 
points on the circle. 

BODY Organizes the flow through VTXCHN. It organizes calls to predict vortex 
shedding, to calculate surface velocities and pressures, and to calculate the 
overall forces and moments on the body. 

CHNGAM Calculates the strengths and locations of vortices shed from a chine cross 
section body. 
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CHNINI 

CMAP 

COMBIN 

CONFOR 

CPPOT 

DFEQKM 

DNUDX 

DPHIDT 

DZDNU 

F 


Calculates the initial strengths and locations of vortices shed from a chine 
cross section body. This is done by approximating the body with chine as 
a circular body with a strake modeled by a vortex lattice. 


Transforms a specified point in the body plane to its corresponding point 
in the circle plane. 

Combines vortices which are separated by a distance less than or equal to 
"RGAM". (The resulting vortex strength is the sum of the individual 
strengths, and its location is the centroid of the combined vortices). 


Calculates the coefficients of the numerical conformal mapping function. 
Transforms the region outside of a polygonal shape with a vertical plane of 
symmetry to the region outside of a circle. 

Calculates a potential pressure distribution at a given x station (the effects 
of separation vortices are not present in the C p calculation). This is used 
in determining a turbulent separation point for the transition region. 

Obtains vortex paths between stations x and x + Ax using a Kutta-Merson 
integration scheme. 

Calculates the noncircular mapping derivative dv/dx. 

Calculates the two-dimensional unsteady pressure term for use in the 
unsteady Bernoulli equation. 

Calculates the differential of the transform of a polygonal shape to a circle. 

Calculates the derivatives of the vortex equations of motion which are 
used in conjunction with DFEQKM to obtain the vortex paths. 
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FDPSDR 

FORCEP 

FPSI 

FPVEL 

INPT 

INVERS 

PLOT 

RDRDX 

RDRDXI 

SEPRAT 

SHAPE 

SKIN 


Calculates the axial velocity due to a three-dimensional source 
distribution. 

Calculates the forces and moments on the body by an integration of the 
surface pressure distribution. 

Calculates the axisymmetric body shape represented by a three- 
dimensional source distribution. 

Calculates the velocity at specified field points. 

Reads in and prints out the body geometry, the flow conditions, various 
program options, and restart quantities. 

Solves simultaneous linear equations, [A]x = b. 

Creates a line printer plot of the body and its vortex wake at any given x 
station. It is composed of the following subroutines: PLOTV, PLOTA2, 
PLOTA5, PLOTA6, PLOTA7, and PLOTA8. 

Calculates r and dr/dx at a given x station. 

Initializes the calculation of r and dr/dx for subroutine RDRDX. 

Predicts the position in the crossflow plane where separation occurs based 
on the Stratford criteria for laminar and turbulent flows. 

Determines the body geometric parameters from a table look-up of input 
data. 


Estimates the axial skin friction force on the body. 
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SORDIS 

SORVEL 

SRCVEL 

SUM 

TRANSN 

UMAP 

VCENTR 

VCNTR 

VCNTRP 

VINFLU 

VLOCTY 


Calculates the source distribution representing an axisymmetric body of 
revolution given the body shape and slope at each x station. 


Calculates perturbation velocities due to a three-dimensional source 
distribution. Used by subroutine SORDIS in calculation of the source 
distribution. 

Calculates perturbation velocities due to a three-dimensional source 
distribution. 

Analytical mapping of an elliptical cross section. 

Calculates the body separation points in the transition region. 

Calculates the axial velocity, 34>/ 3x, due to noncircular mapping and two- 
dimensional singularities. 

Calculates the centroid of each vortex field. 

Calculates the centroid of any set of vortices. This routine is used for the 
d^/dt calculation. 

Calculates the actual centroid of vorticity for subroutines VCNTR and 
VCENTR. 

Calculates the influence of a unit strength vortex and its image. 

Calculates the velocity components in the cross flow plane at a specified 
field point. 
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VIABLE Sets up the arrays of free vorticity for two cases: (1) The vorticitv in the 
working arrav CAM is distributed into the individual arrays GAMP, 
CAMM, GAMR, GAMA, or (2) the vorticity in the individual arrays GAMP, 
GAMM, GAMR, GAMA is distributed into the working array, GAM, in the 
proper order. 

Z Calculates the transformation of a polygonal shape to a circle of radius RZ. 

Program Limitations and Suggestions 

Program VTXCHN predicts the aerodynamic forces acting on slender bodies alone at 
high angles of attack, with some limitations. Certain of the limitations that follow are 
basically rules-of-thumb for the use of the program which have been determined through 
program verification and experience with the code. 

Ghinp cross section bodies. - There are several things that are required when modeling 
bodies with chine cross sections. First, because of limitations in the numerical conformal 
mapping procedure, the input geometry specification of the chine cross section may have to 
be simplified. If the numerical mapping fails to converge or if the message "ERROR IN THE 
SUM OF EXTERIOR ANGLES" occurs, the number of points specifying the cross section 
may need to be reduced while maintaining the general shape of the cross section. 

The source distribution input for a chine body must correspond to an equivalent area 
distribution body. The equivalent area distribution radius, r eq/ must be calculated external 
to the program. The mapped radius, r Q , is also required by the program. The only method 
to obtain the mapped radius is to run the program. It is only necessary to run the code for 
the first axial station (i.e. set XI = XF) since all of the mapping information is calculated at 
the beginning. The input requires that the mapped radius be input with the same resolution 
as the equivalent radius; therefore, this may require an effort or utility program external to 
VTXCHN to obtain r Q at the axial station specified for r^. 
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With a chine cross section it is possible for tracking problems to occur, and a STOP 40 
will occur. The chine cross sections are more prone to this problem than circular or elliptic 
cross sections. If this problem is encountered, the following variables can be altered to 
alleviate the tracking problems: 

E5 tolerance for vortex tracking 

RCORE vortex core size normalized by body diameter 

DX axial spacing between tracking stations 

Mach number .- The major limitation in the program is the Mach number or 
compressibility correction. The compressibility correction added to VTXCHN is based on 
small perturbation theory which restricts the body to be slender, the angle of attack to be 
small, and the Mach number to be approximately 0.8 or less. Program VTXCHN including 
the compressibility correction models subsonic flow with linear theory and is not capable of 
modeling mixed subsonic/supersonic flows or shock waves. These restrictions result from 
the assumptions made for small perturbation theory, and any cases outside these 
restrictions should be viewed with caution; however, VTXCHN has been applied to blunt 
nose shapes such as spheroids or ellipsoids with good success. Typical missile shapes with 
slender noses pose no problems to the method. 

Source di stribution - The following suggestions should be observed when using the 
source distribution calculation in VTXCHN. The source distribution in VTXCHN places 
three-dimensional point sources on the axis of the body (or equivalent body in a noncircular 
case). For subsonic Mach numbers, it is necessary to close the body to calculate a smooth 
source distribution. Generally, the body can be closed at its base by including a smooth 
boattail, and a successful technique for most cases is to close the model with its nose. The 
source distribution does not handle rapid changes or discontinuities in body slope well. 

The source distribution calculation often has a problem with short pointed noses as 
well. Short blunted noses seem to cause no difficulty because they resemble the nose of a 
Rankine body, a single point source in a uniform stream. Sharp pointed noses often cause 
an oscillation in the source strength on the afterbody resulting in an oscillation in the 


- 52 - 



velocities and thus an oscillation in the axial pressure distribution. This has been observed 
with a 2-caliber ogive nost but was not apparent for a 2.5-caliber ogive. The oscillations 
were minor in their effect on the total calculation and result, but the user should be aware of 
their possibility. 

Incidence angle .- The well-known and little-understood phenomenon of wake 
asymmetry at large incidence angles (or c > 30 ) is available as an option in VTXCHN. As 
described in Reference 1, the naturally symmetric nature of the vortex shedding analysis 
must be perturbed to produce an asymmetric solution. The user must be aware of the 
limitations associated with this portion of the code in order to apply appropriate judgement 
in the evaluations of the results. Based on the author's experience with asymmetry, care 
must be taken in the choice of the perturbation parameter and its effect on the results. This 
is to serve as a warning, and the user should be aware that this portion of the code is in need 
of additional development. 

Transition .- The transition capability included in the code for smooth bodies is in an 
early stage of development. The character of the vortex wake, whether originating from 
laminar or turbulent separation points, has an important effect on the predicted nonlinear 
aerodynamics characteristics. As previously discussed (Ref. 1), additional development is 
needed to better understand the origin and length of the transition regions. The user is 
cautioned to be very careful in the interpretation of predicted results in flow regimes which 
may be in a transition region as defined by Lamont in Reference 33. 

The transition option in the program is not applicable to chine cross sections. With 
chine cross sections, separation is fixed at the sharp edges. 

Secondary separation .- The capability to predict secondary separation or separation in 
the reverse flow region on the lee side of the body is included in VTXCHN. As discussed in 
a later section, there are certain flow conditions in which secondary separation has an 
important role; for example, asymmetric shedding at large incidence angles. Numerical 
problems associated with tracking the secondary vorticity have been observed, and some 
additional effort is required to make secondary separation a useful option. Therefore, it is 
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recommended that this option be treated as a research feature and not used for production 
calculations. 


Error Messages and Stops 

The code VTXCHN has numerous internal error messages. These messages are 
generally explicit and are described below. There are several execution terminations at 
numbered stops within the program. These are described in the table below with some 
suggestions to correct the problem. 

The error message "CONVERGENCE NOT OBTAINED IN INTEGRATION" is printed 
out when the trajectory calculation does not converge within the specified tolerance, E5. 
Also printed out is the axial station, XI, and the integration interval, H. This error results in 
a STOP 40 which is discussed in the table to follow. 

In the calculation of the source distribution it is possible to get the message "MATRIX 
IS SINGULAR. This occurs when one of the pivot elements in the influence matrix is zero. 
This error results in a STOP 300. 

The other major error messages are found in the numerical conformal mapping routine 
CONFOR. The message "‘"ERROR IN THE SUM OF EXTERIOR ANGLES*’’" is printed if 
the routine has difficulty calculating appropriate body angles within a specified tolerance. 

Another message "’’’COEFFICIENT AN(1) = GREATER THAN +- 0.0001 AFTER 

ITERATIONS OUT OF " is printed on the rare occasions when the mapping procedure is 

having difficulty converging on a solution. These messages are informational and do not 
result in a program STOP. However, the user is advised to check the transformed shape 
carefully with the actual shape before using the results. 

Termination 

Message Description 

STOP 1 This is the normal stopping condition. Program VTXCHN 
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STOP 40 


STOP 41 


STOP 44 


STOP 45 


STOP 46 


The vortex trajectory integration has failed to converge at an axial 
station. This problem is usually caused by: (1) two or more vortices 
in close proximity rotating about one another so rapidly that the 
integration subroutine cannot converge on a solution, or (2) a vortex 
is too near the body surface. These situations are often dependent 
on specific flow conditions and the body shape. The following steps 
should be tried to resolve the problem. Increase the integration 
tolerance E5, or increase the size of the combination parameter 
RGAM so that it is greater than the distance between the 
troublesome vortices. It often helps to increase the integration 
interval DX and restart the calculation procedure from a previous 
station. Subroutine BODY 

The vortex trajectory calculation has resulted in a vortex inside the 
body This is a rare but serious problem which indicates a major 
problem with the solution. Subroutine BODY 

The total number of shed vortices exceeds 200. To correct this, 
combine some vortices and restart the calculation, or increase the 
integration interval, DX, so that fewer vortices are shed. Subroutine 
BODY 

The maximum number of positive or negative primary separation 
vortices exceeds 70. Use the same remedy as for STOP 44. 
Subroutine BODY 

The maximum number of positive or negative secondary (reverse 
flow) separation vortices exceeds 30. Use the same remedy as for 
STOP 44 but only combine like vortices. If one of the above steps 
does not work, the reverse flow calculation can be turned off by 
setting NSEPR = 0. Subroutine BODY 
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STOP 201,202 1 here is an error in the vortex tables in subroutine VTABLE. Check 

the input values of NBLSEP, NSEPR, NVP, NVM, NVA, NVR, and 
RGAM in Items 1 through 6. 

STOP 300 There is a singular influence matrix in the solution for the three- 

dimensional source distribution. Check the input geometry for 
errors. Subroutine INVERS 

STOP 301 The maximum number of iterations is exceeded in the source 

calculation. Check the input geometry. This problem can occur if 
the body has a long slender nose followed by a cylindrical section 
and then a long slender tail. If this is the case, one solution is to 
close the body off with a blunter and/or shorter tail. Subroutine 
SORDIS 


Input Description 

This section describes the input data and format required by program VTXCHN. The 
input formats are shown in Figure 22, and the definitions of the individual variables are 
provided in the following paragraphs. Note that some length parameters in the input list 
are dimensional variables; therefore, special care must be taken that all lengths and areas are 
input in a consistent set of units. 

The remainder of this section includes a description of each of the input variables and 
indices required for the use of program VTXCHN. The order follows that shown in Figure 
22 where the input format for each item and the location of each variable on each input line 
is presented. Data input in I-format are right justified in the fields, and data input in F- 
format may be placed anywhere in the field and must include a decimal point. 

Item 1 is a single input line containing 15 integer flags. These flags are right justified in 
a five column field and are used to specify several options available in VTXCHN. 
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I ( 1515 ) NC1R Cross section flag, determines the type of body cross 

section. 

= 0, circular body 

= 1, elliptical body 

= 2, arbitrary body, similar cross section at all axial 
stations 

= 3, arbitrary body, nonsimilar cross sections 


NCF Mapping flag, specifies whether the required mapping 

is input or calculated. (NCIR = 2 or 3) 

= 0, calculate mapping 

= 2, input mapping 


ISYM Vortex shedding symmetry flag. 

= 0, symmetric shedding, (assumes symmetry to save 
computation time) 

= 1, nonsymmetric shedding, (if PHI * 0 in item 5 
ISYM = 1) 


NBLSEP Separation flag. 

= 0, no separation 

= 1, laminar separation 

= 2, turbulent separation 

NSEPR Reverse flow or secondary separation flag. 

= 0, no reverse flow separation (preferred) 
= 1, reverse flow separation 


NDFUS Vortex core model flag. 

= 0, potential vortex model 
= 1, diffuse vortex model (preferred) 
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NDPHI 


1NP 


NXFV 


NFV 


NVP 


NVR 


NVM 


NVA 


NASYM 


I vvo-dimensional unsteady pressure term flag. 

= 0, omit 30 / dx term in C p equation 
= 1, include 30/ dx term (preferred) 

Nose force flag. 

= 0, estimate nose forces ahead of the first axial 
station (preferred) 

= 1/ zero nose forces ahead of the first axial station 
(for restart calculation) 

Number of x stations at which the flow field is 
calculated and/or special output is generated. 
(0 £ NXFV £ 8) 


Number of field points for flow field calculation. 
(0 £ NFV ^ 500) 


Number of +T vortices on +y side of body for a restart 
calculation. (0 £ NVP £ 70) 

Number of reverse flow vortices on +y side of body for 
a restart calculation. (0 £ NVR £ 30) 

Number of -r vortices on -y side of body for a restart 
calculation. (0 £ NVM £ 70) 

Number of reverse flow vortices on -y side of body for a 
restart calculation. (0 £ NVA £ 30) 

Asymmetric vortex shedding flag, (not used on chine 
portion of body) 

= 0, no asymmetry 

= 1, forced asymmetry (not recommended) 
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Item 2 is a single input line containing 3 integer flags. I he first flag NVPHI controls 
how d<p/dx is calculated for the vortices, the second flag specifies whether the body has a 
chine shaped section, and the third flag specifies whether the dA k / dx term of Equation (9b) 
is included in dv/dx. 


Item Variable Descript ion 

2 (315) NVPHI d(p/dx flag, specifies how vortices are treated in the 

d<f>/dx calculation. 

=-l, all vortices are considered independent for the 
d<p/dx calculation (preferred) 

> 0, the most recent NVPHI vortices (on each side of 
the body) are considered independent for the 
d<t>/dx calculation, and the remaining vortices 
are included by a centroid representation. The 
NVPHI vortices considered independently can 
be considered as the feeding sheet. 

NVPHI =-l is suggested. 


NCHINE chine flag, determines whether the body has a chine 

shaped cross section. 

= 0, no chine section, body determined by NCIR 
= 1, chine section for all or part of the body 


NDADT flag which determines whether the term in 

Equation (9b) involving dA k /dx is included in the 
calculation of dv/dx. dA k /dx is calculated numerically 
by the program. 

= 0, dA k /dx term not included in dv/dx 
= 1, dA k /dx term is included in dv/dx 
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Item 3 is a single input line containing 10 integer flags/ parameters. These values are 
used to specify several options available in VTXCHN. 


1 tern Variable 


Description 


3(1015) NHEAD Number of title lines in Item 3. 

(NHEAD > 0) 


NPRNTP Pressure distribution print flag. 

= 0, pressure distribution output only at special 
output stations (NXFV > 0) 

= 1, pressure distribution output at each x station 

NPRNTS Vortex separation print flag. 

= 0, no separation output 

= 1, separation point summary (preferred) 

= 2, detailed separation calculation (for diagnostic 
purposes only) 


NPRNTV V ortex field summary output flag. 

= 0, vortex field is printed at special output stations 
only (NXFV > 0) 

= 1, vortex field is printed at every x station 

NPLOTV Vortex field printer-plot flag. 

= 0, no plots 

= 1, plot full cross section, constant scale 

= 2, plot half cross section, constant scale 

= 3, plot full cross section, variable scale 
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N PLOT A 


Additional output t lag. Option to have pressure 
distribution, vortex lield summary, and a plot of the 
vortex field regardless of the above set options. 

= 0, no additional output 

= 1, additional output only at special x stations 
(NXFV > 0) 

= 2, additional output at each station 

NTH Numerical conformal mapping angle flag. 

= 0, calculate circle 0's for mapping at 5" increments 
(preferred) 

= 1, input d's 

NCORE Vortex core radius flag. 

= 0, maximum allowable local vortex core size is 

.05*D 

= 1, no upper limit on the local vortex core size, the 
vortex core size is given by RCORE*D 

NSKIN Axial skin friction drag estimation flag. 

= 0, no skin friction estimation 
= 1, axial skin friction estimated 

NCOMP Compressibility correction flag. 

= 0, incompressible calculation 

= i, compressibility correction applied (Mach 

number must be input in Item 6) 

item 4 is a set of NHEAD title/summary lines which are printed at the top of the first 
output page. They are reproduced in the output just as they are input. 
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Item 

Variable 

Description 

4 (20A4) 

TITLES 

NHEAD title/summary lines to be output on first page. 

Item 5 is a single input 
aerodynamic coefficients. 

line containing reference information needed to form the 

Item 

Variable 

Description 

5 (5F10.5) 

REFS 

Reference area, S ref . REFS > 0. 


REFL 

Reference length,l ref . REFL > 0. 


XM 

Axial position of the moment center, measured from the 
body nose. 


SL 

Body length, L. 


SD 

Maximum body diameter, D. 

Item 6 is 

a single input line containing the flow conditions. 

Item 

Variable 

Description 

6 (3F10.5) 

ALPHAC 

Angle of incidence (degrees). (0‘Sa c < 90°) 


PHI 

Angle of roll, <f> (degrees). For chine bodies the user 
should be aware that separation is forced at the sharp 
edge of the chine regardless of the value of PHI. It is the 
user's responsibility to determine the limit of PHI for a 
given chine shape. 
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RE Reynolds number based on maximum body diameter, 

(V.D/V). 

MACH Free-stream Mach number. 


Item 7 is a single input line which specifies the axial extent of the run, the transition 
region, and parameters associated with the vortex wake 


Item 

Variabe 

Description 

7 (8F10.5) 

XI 

Initial x station, XI > 0, dimensional quantity. (Note: 
XI 2: DX/2) 


XF 

Final x station, XF > XI, dimensional quantity. 


DX 

x increment for vortex shedding calculation. Typical 
value, DX * D/2. For bodies with chine cross sections it 
is recommended that DX be between 0.10*D and 0.25 D 


XTRl 

Beginning of transition region, dimensional quantity. 
This option is not available on a chine portion of the 
body. 


XTR2 

End of transition region, dimensional quantity. This 
option is not available on a chine portion of the body. 
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hMKF Minimum shed vortex starting distance from the body- 

surface. This value is not used on a chine portion of the 
body. 

= TO, shed vortices located according to stagnation 
point criterion 

> 1.0, minimum radii position of shed vortices. 
Typical value, EMKF = 1.05 

RGAM Vortex combination factor. 

= 0.0, vortices are not combined (preferred) 

>0.0, radial distance within which vortices are 
combined 


Vortex reduction factor, 6, to account for observed 
decreases in vortex strength. 

= 0.6, for bodies with bases 

= 1 .0, for closed bodies 

If NCHINE = 1, this parameter is not used. 


Item 8 is a single input line containing the integration tolerance, asymmetry 
perturbation information, alternate separation criteria, and the core radius. 


Item Variable 


Description 


8 (8F10.5) E5 


Error tolerance for the vortex trajectory calculation. 
Typical range, E5 = 0.01 to 0.05, or larger if appropriate. 


XTABL 


Location at which a table of points on the body and on 
the circle is printed. 

= 0.0, no table of points is printed (preferred) 

> 0.0, table of points for all x < XTABL 
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XASYIV I 


Initial x station for asymmetric perturbation. This 
option is not available on a chine portion of the body. 

XASYMF Final x station for asymmetric perturbation. This option 

is not available on a chine portion of the body. 

DBETA Angular displacement of separation points for 

asymmetry perturbation. This option is not available on 
a chine portion of the body. 

Sppp Stratford laminar criterion, (not used on chine portion of 

body) 

= 0.0, program uses SEPL = .087 (preferred) 

> 0.0, program uses input value 

SEPT Stratford turbulent criterion, (not used on chine portion 

of body) 

= 0.0, program uses SEPT = .350 (preferred) 

> 0.0, program uses input value 

RCORE Ratio of the local core radius to the local body diameter. 

Default value is .025. Maximum allowable value is 
,025*D/r o . 

Item 9 is a single input line containing two flags for the source distribution. The source 
distribution may be input (NSOR > 0), calculated from the body geometry (NSOR = -1), or 
omitted from the solution (NSOR = 0). If the Mach number is not equal to zero (NCOMP > 
0), then a different source distribution is required for each Mach number. 
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Description 

Source distribution flag. 

=-1, calculate the source distribution from the 
geometry input 

= 0, no sources (not recommended) 

> 0, number of sources to be input 

Source distribution print flag. 

= 0, output source locations and strengths 
= 1, above output plus input and source geometry 
comparison 

= 2 , above output plus source-induced surface 
velocities 

Omit Items 10 and 1 1 if NSOR s 0. 


Items 10 and 11 contain the nondimensional source locations and nondimensional 
source strengths, respectively. 


Item 

Variable 

Description 

10 (6E12.5) 

XSRC 

Nondimensional source locations, x/L. (NSOR values, 6 
per input line) x/L may be greater than 1 .0. 

Item 

Variable 

Description 

11 (6E12.5) 

QS 

Nondimensional source strengths, Q/4jrL 2 V oo cosa c . 

Items 12 through 15 contain the body geometry information. This information must be 

included for all bodies whether circular or noncircular. 


Item 

Variable 

9 (215) 

NSOR 


NPRT 
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Item 


Variable 


Description 


12(15) NXR 


Number of entries in the body geometry table. 
(1 s NXR ^ 101) 


Item Variable 


Description 


13 (8F10.5) XR 


Nondimensional axial stations in the geometry table, 
x/L. (NXR values, 8 per input line) 


Item Variabl e 


Description 


14 (8F10.5) R 


Nondimensional equivalent area body radius at the 
NXR x/L-stations, r^/L (NXR values, 8 per input line). 


Item Variabl e 


Descriptio n 


15 (8F10.5) DR 


Body slope at the NXR x/L-stations, dr eq / dx. (NXR 
values, 8 per input line) 


Omit Items 16 and 17 if NC1R * 1- 


Items 16 and 17 contain the nondimensional horizontal and vertical half axis lengths if 
the cross sections are elliptical, NCIR = 1. 


Item Variab le Description 


16(a) (8F10.5) AE 


Nondimensional horizontal half axis length of the 
elliptical cross section at the NXR x/L-stations, a/L. 
(NXR values, 8 per input line) 
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Item 


Variable 


Description 


16(b) (8F10.5) BE 


Nondimensional vertical half axis length of the elliptical 
cross section at the NXR x/L-stations, b/L. (NXR 
values, 8 per input line) 


Item Variable 


Description 


17(a) (8F10.5) DAE 


Slope of horizontal meridian of the elliptical cross 
section at the NXR x/L-stations, da/dx. (NXR values, 8 
per input line) 


Item Variable 


Description 


17(b) (8F10.5) DBE 


Slope of vertical meridian of the elliptical cross section 
at the NXR x/L-stations, db/dx. (NXR values, 8 per 
input line) 


Omit Items 18(a) and 18(b) if NCIR = 0. 


Items 18(a) and 18(b) contain the mapped body geometry information at the NXR x/L- 
stations specified in Item 13. 


Item Variable 


Description 


18(a) (8F10.5) RO 


Nondimensional transformed (mapped) body radius at 


the NXR x/L stations, r Q /L (NXR values, 8 per input 
line). For an ellipse, r Q = (a + b) /2. 
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Item 


Variable 


Description 


18(b) (8F10.5) DRO Body slope at the NXR x/L-stations, dr Q / dx. (NXR 

values, 8 per input line) 

Omit Items 19 through 27 i: NC1R < 2. 

Items 19 through 27 are used to specify a body of arbitrary cross section shape. This 
body should also be polygonal in shape (i.e., negative interior angle can cause problems in 
the numerical conformal mapping procedure). 

Omit Item 19 if NTH = 0. 

Item Variable ? Description 

19 (8F10.5) TH Special table of circle angles for arbitrary body mapping, 

(i = 1, 73), to be used only if the density of points on the 
cross section needs adjusting to improve the resolution 
of pressure distribution. Not recommended under 
normal conditions. 


Item Variable Description 

20(a) (215) MNFC Number of mapping coefficients for the arbitrary body. 

(MNFC ^ 100). Generally good mappings can be 
obtained with 20 to 30 coefficients. 


MXFC Number of axial stations at which the arbitrary body is 

specified. (1 £ MXFC ^ 10). For a similar shape at all 
cross sections, MXFC=1. 
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Item 


Variable 


Description 


20(b) (8F10.5) XFC x stations at which the arbitrary body is specified. 

(MXFC values, 8 per input line). For a similar cross 
section body, MXFC=1, XFC(l) i XI. 

Omit Item 21 if NCHINE s 0 

Item 21 is a set of MXFC integers to indicate which input corner point of each input 
cross section is to be considered as the chine edge. Figure 23(b) illustrates a chine and the 
choice of NCSEP. The input corner points are read in Items 22 and 23. (MXFC values, 10 
per input line). 

Item Variable Description 

21 (15) NCSEP Number of the input corner point defining the chine 

edge. Corner points are read in items 22 and 23. It is 
recommended that NCSEP be the same for all input 
chine cross sections. 0 s NCSEP ^ NR. Omit item 21 if 
NCHINE s 0. 


Omit Items 22 and 23 if NCF > 0. 

Items 22 and 23 are repeated for each of the MXFC axial stations when the 
transformation coefficients are to be calculated (NFC = 0). 
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Item 

Variable 

Description 

22 (15) 

MR 

Number of points describing the arbitrary body cross 
section (2 < NR ^ 30). (Note: NR must be the same for 
all x stations.) 

Item 

Variable 

Description 

23(a) (6E12.5) 

XRC 

Horizontal coordinates of the arbitrary body cross 
section. (NR values, 6 per input line). The convention 
for ordering the coordinates from bottom to top in a 
counter-clockwise fashion is observed. Figure 23. This 
corresponds to inputting the +y side geometry. 
Right/ left body symmetry is required. 

23(b) (6E12.5) 

YRC 

Vertical coordinates of the arbitrary body cross section. 
(MR values, 6 per input line). See Figure 23. 


Omit Items 24 through 27 if NCF = 0. 

Items 24 through 27 are repeated as a group for each of the MXFC axial stations if the 
numerical transformations are available from a previous calculation. The only purpose for 
providing the option to input the mapping parameters is to eliminate the need to recompute 
the numerical mapping in subsequent runs and save computation time. 


Item 

Variable 

Description 

24 (2F10.5) 

zzc 

Mapping offset distance. 


RZC 

Radius of transformed cross section. 
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Item 


Variable 


Description 


25 (15) 


NR 


Number of points describing the arbitrary body cross 
section (NR ^ 30). (Note: NR must be the same for all x 
stations.) 


Item 


Variable 


Description 


26(a) (6E12.5) XRC 


Horizontal coordinates of the arbitrary body cross 
section. (NR values, 6 per input line). See Figure 23 and 
convention for Item 23. 


26(b) (6E12.5) YRC 


Vertical coordinates of the arbitrary cross section. (NR 
values, 6 per input line). See Figure 23. 


26(c) (6E12.5) THC 


Circle angles, 9 , of body points. (NR values, 6 per input 
line) 


Item 


Variable 


Description 


27 (6E12.5) AFC 


Mapping coefficients. (MNFC values, 6 per input line) 


The derivative of the mapping required in Equation (9b) is calculated by the program. 


Omit Items 28 - 29 if NCHINE < 0. 


Item 28 is a single input line containing information required to start the solution by 
approximating the chine as a thin strake on a circular body. 
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Item 


Variable 


Description 


28 (I5,F10.5) NSTRTC number of axial stations at which the vortex lattice 

approximation of the chine is used to start the solution. 
A value between 0 and 5 generally works the best, but 
the user is encouraged to investigate the effect of this 
parameter. 

FVCHN fraction of the exposed strake span where the boundary 

condition is applied to determine the initial chine 
strength using a vortex lattice approximation. 
Suggested value 0.50. 

Item 29(a) is a single input line containing the' starting and ending axial stations of the 
chine portion of the body. 

Item Variable Description 

29(a) (2F10.5) XCHNB axial station at which chine separation begins. 

XCHNE axial station at which chine separation ends. 

Item 29(b) is a single input line containing three values of NCIR for a chine body. See 
input Item 1 for a description of NCIR. 

Item Variable Description 

29(b) (315) NCH(l) value of NCIR for portion of body in front of axial 

station XCHNB. Only NCIR = 0, 2 ,or 3 is available here. 
If the body has an elliptical nose section, it should run as 
a noncircular body (NCIR = 2 or 3). 
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NCH(2) value of NCIR for portion of body between axial station 

XCHNB and XCHNE. For the chine section, NCIR is 2 
or 3. If NCH(l) = 0 and NCH(2) = 3, then the first input 
body cross section should be a circle with a radius equal 
to the body radius at the start of the chine section. This 
is required to ensure transition from the circluar 
portion of the body to the chine portion. 

NCH(3) value of NCIR for portion of body after axial station 

XCHNE. Only NCIR = 0, 2 ,or 3 is available here. If the 
body has an elliptical afterbody section, it should run as 
a noncircular body (NCIR = 2 or 3). 

Item 29(c) is a single input line containing two values of NBLSEP for a chine body. See 
input Item 1 for a description of NBLSEP. 

Item Variable Description 

29(c) (215) NCBL(l) value of NBLSEP for portion of body in front of axial 

station XCHNB. 

NCBL(3) value of NBLSEP for portion of body after axial station 

XCHNE. 

The value of NCBL(2) is set in the code but is not used. 

Item 29(d) is a single input line containing two values of VRF for a chine body. See 
input Item 7 for a description of VRF. 
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Item 


Variable 


Description 


29(d)(2F10.5) VRFC(l) value of VRF for portion of body in front of axial station 

XCHNB. See input Item 7 for recommended values. 

VRFC(3) value of VRF for portion of body after axial station 

XCHNE. See input Item 7 for recommended values. 

The value of VRFC(2) is set in the code but is not used. 


Omit Items 30 and 31 if NXFV = 0. 

Item 30 is a single input line containing the x stations at which additional output is 
printed and/or calculated if NXFV > 0. 


Item Variable 


Description 


30 (8F10.5) XFV 

Omit Item 31 if NFV = 0. 


x stations at which additional output is printed or 
calculated. (NXFV values, 8 values maximum) 


Item 31 is a set of NFV input lines which contains the coordinates of the field points 
where the velocity components are to be calculated at each of the axial stations in Item 30. 


Item Variable Description 

31 (2F10.5) YFV,ZFV y and z coordinates of the field points for the velocity 

field calculation. (NFV input lines, maximum of 500). It 
is recommended that all points lie outside the body 
surface at all axial stations. 
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Omit Items 32 through 37 if NVP + NVR + NVM + NVA = 0 


Items 32 through 36 provide information needed to restart the calculation procedure. It 
is not necessary to input all four types of vorticity in a restart calculation. Note that 
program VTXCHN write a file to FORTRAN Unit 1 which contains the required input for 
Items 32 through 36. 

Item 32 is a single input line containing the nose force and moment coefficients at the 
restart point; that is, X = XI. These values may be set equal to zero, but the resulting forces 
and moments from the current calculation will be in error. 


Item 

Variable 

Description 

32 (5F10.5) 

CN 

Normal- force coefficient 


CY 

Side-force coefficient 


CA 

Axial pressure force coefficient 


CDI 

Axial skin friction force coefficient 


CM 

Pitching-moment coefficient 


CR 

Yawing-moment coefficient 


CSL 

Rolling-moment coefficient 


Omit Items 33 if NVP = 0. 


Item 33 is a set of NVP input lines which specify a cloud of positive primary body 
separation vorticity, usually on the right side of the body. The variable XSHEDP associated 
with each vortex is used only to identify individual vortices and permit the user to follow 
individual vortices during the calculation. 
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Item 

Variable 

Description 

33 (4F10.5) 

GAMP 

+y primary vorticity, T/V^. 


YGP 

Horizontal coordinate, y, of vortex. 


ZGP 

Vertical coordinate, z, of vortex. 


XSHEDP 

Origin of the vortex. 


Omit Item 34 if NVR = 0. 

Item 34 is a set of NVR input lines which specify a cloud of reverse flow (secondary) 
body separation vorticitv. This array is a convenient way to put an arbitrary cloud of 
additional vorticity which is to be maintained separately from the other body vorticity in 
the field. 


Item 

Variable 

Description 

34 (4F10.5) 

GAMR 

+y reverse flow vorticity, r/V^. 


YGR 

Horizontal coordinate, y, of vortex. 


ZGR 

Vertical coordinate, z, of vortex. 


XSHEDR 

Origin of the vortex. 


Omit Item 35 if ISYM = 0 or NVM = 0. 

Item 35 is a set of NVM input lines which specifies a cloud of negative primary body 
separation vorticity, usually on the left side of the body. In the case of a symmetric flow 
field (ISYM=0), this block of vorticity is automatically defined as the mirror image of the 
positive vorticity input in Item 33. 
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Item 

Variable 

Description 

35 (4F10.5) 

GAMM 

-y primary vorticity, r/V . 


YGM 

Horizontal coordinate, y, of vortex. 


ZGR 

V ertical coordinate, z, of vortex. 


XSHEDR 

Origin of the vortex. 

Omit Item 36 if ISYM = 0 or NVA 

= 0. 

Item 36 is a set of NVA input lines which specifies a cloud of reverse flow (secondary) 
body separation vorticity. This block of vorticity is analogous to Item 34, and it is omitted if 

NVA = 0 or if ISYM = 0. As with 
defined as the mirror image of 
required by ISYM = 0. 

the previous item, this block of vorticity is automatically 
the negative vorticity input as Item 34 if symmetry is 

Item 

Variable 

Description 

36 (4F10.5) 

GAMA 

-y reverse flow vorticity, r/V . 


YGA 

Horizontal coordinate, y, of vortex. 


ZGA 

Vertical coordinate, z, of vortex. 


ZSHEDA 

Origin of the vortex. 


Item 37 contains the pressure distribution at the initial x station for a restart calculation. 
These values should be obtained from a previous run. If a symmetric case is being 
considered, ISYM=0, only half of the pressure distribution is needed (37 points), for a 
asymmetric case 73 points are required. 


Item 

Variable 

Description 

37 (3F10.5) 

YBO 

y body coordinate. 


ZBO 

z body coordinate. 


CPO 

pressure coefficients at (YBO, ZBO). 
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This concludes the input deck for a single run of VTXCHN. The code is not set up to 
stack input for multiple cases because of the sometimes long and often unpredictable run 
times. 


Input Preparation 

In this section, some additional information is provided to assist the user in the 
preparation of input for certain selected problems. The previous section on the input 
description must be used to understand the individual variables which go into VTXCHN, 
and this section will permit the user to select the appropriate input to get optimum results 
from the code. 

Numerical mapping .- The specification of the appropriate numerical mapping 
parameters (Items 18 through 28) depends on the shape of the cross section. In the interest 
of optimum computation time, the fewest possible terms in the transformation series 
(MNFC in Item 19) should be used which will produce an adequate mapping. Each 
different shape should be checked by the user to determine the proper number of 
coefficients. For example, a simple cross section which is similar to a circle may require as 
few as 10 coefficients. Others may require as many as 30 coefficients for an excellent 
mapping. If a large number of calculations are to be made for a noncircular shape, it is 
worthwhile for the user to try different numbers of coefficients in an effort to find an 
optimum number. Not only should the shapes be compared carefully, the unseparated flow 
pressure coefficients should be compared to evaluate the quality of the mapping. 

The numerical mapping now included in VTXCHN has been used for a variety of cross 
section shapes, and some problems have been observed. One general problem area is bodies 
with concave sides. The mapping procedure has difficulty converging on a set of 
coefficients. Usually the body shape can be modified slightly, to give the concave region a 
flat or slightly convex shape, and the mapping converges easily. This small change in the 
shape probably has very little effect on the predicted pressure distribution. This is an 
example of the type of problem for which the user should be alert. For chine cross sections. 
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several problems can occur due to concave or drooped sides. The user may need to be 
ceative in coming up with the best approximation of the real cross section. This may require 

specifying fewer points to describe the body cross section in order to obtain a good 
numerical mapping. 


I ntegration interval .- Choice of the appropriate integration interval (DX in Item 7) 
depends to a large extent on the type of results of interest. If the user is interested in 
detailed nose loads and separation on the nose where the radius is changing rapidly, the 
interval should be smaller than the recommended D/2. On very long bodies, the 
recommended interval is usually adequate for accurate predictions. For chine bodies the 
interval is recommended to be smaller than D/4 and as small as D/10, but smaller intervals 
can cause numerical problems. 


V ortex core .- The core size specification is referenced to the local body diameter 
(RCORE in Item 8). The default value is 0.025, and this value is adequate for most cases. 
The larger the core radius specified, the more the vortex effects are reduced. To eliminate 
the possibility of negating the vortex effect entirely, a maximum vortex core of .05 times the 
maximum body diameter is built into the code. Except in unusual situations, the default 
value is recommended. 


Chine Sample Cases 


In this section, two sample input cases are described to illustrate the use of program 
VTXCHN to model bodies with chine cross section shapes. Nine sample case inputs dealing 
with circular, elliptic, and noncircular bodies are presented in Reference 1. The purpose of 
the following sample cases is to point out the available features and options for chine 
shaped bodies, and to provide a range of sample inputs which will help users prepare input 
for specific cases. 

Sample Case ! _■- Chine Sample Case 1, Figure 24, is a test run based on the chine body 
in Reference 32. This run will also be used as the sample output case in the next section. It 
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is recommended that this case be run initially to provide a check with the results presented 
in a following section. 

Input Item 1 indicates that the case consists of a body with noncircular/ nonsimilar 
cross sections, NCIR = 3, and a symmetric vortex wake, ISYM = 0. The mapping flag 
NCF = 0 indicates that the numerical mapping is calculated. The separation criteria is 
laminar, NBLSEP = 1, but no reverse flow separation is to be computed, NSEPR = 0. A 
diffuse vortex core model is to be used, NDFUS = 1, and the two-dimensional unsteady 
pressure term is included in pressure calculation, NDPHI = 1. INP = 0 indicates that the 
nose forces ahead of the first axial station will be estimated. NXFV = 0 and NFV = 0 indicate 
that no special output stations or field points are to be read in Items 30 and 31. No inital 
vortices are to be read, NVP = NVR = NVM = NVA = 0, and asymmetric shedding is not 
forced, NASYM - 0. 

Input Item 2 contains three additional flags. NVPHI = -1 indicates that all vortices are 
considered independent for the d<*/dt calculation, NCHINE = 1 indicates that the body has 
a chine cross section portion, and NDADT = 0 indicates that the dA k /dx term in 
Equation (9b) is omitted. The chine flag will require further input which will overide some 
of the program control flags in input item 1. When a body with a chine section is being 
investigated, the forced asymmetry (NASYM) option is not available. For chine bodies the 
values of NCIR and NBLSEP are reset by information read in input Items 29(b) and (c). 

Input Item 3 indicates there will be four title cards input in Item 4. NPRNTP — 1 
specifies the pressure distribution be output at all x stations, and NPRNTS = 1 causes a 
summary of the separation points to be output at all x stations. The flag NPRNTV — 1 
results in a summary of the vortex field at all axial stations. NPLOTV = 3 and NPLOTA — 2 
specifies that a line printer plot of the vortex field is created at all x stations if separation has 
occurred. NTH = 0 indicates that the circle 0's are calculated every 5 degrees by the 
program, and NCORE = 0 set the maximum core size to 0.05*D. NSKIN = 0 suppresses 
estimation of the skin friction on the body, and NCOMP = 0 specifies an incompressible 
calculation. 
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Input Item 4 is a set of four title/summary cards which are printed at the top of the 
first output page. 

Input Item 5 contains the reference information for the run. The reference area is 
REFS = 22.145, the reference length is REFL = 8.70, and the moment center, XM, is located at 

station 23.00. The length of the body, SL, is 73.40, and the maximum body diameter, SD, is 
8.70. 

Item 6 contains the flow conditions. The angle of attack, ALPHAC, is 20 degrees, and 

there is no roll angle, PHI = 0.0. The Reynolds number, RE, is 1000000, and the Mach 
number is 0.0. 

Item 7 specifies run parameters. The calculation procedure begins at x station, XI, 0.50 
and continues to station, XF, 14.5 in DX increments of 0.50. A transition region is not 
specified, XTR1 = XTR2 = 0.0. Separation vortices are initially placed a minimum 1.05 radii, 
EMKF, from the axis of the missile. Vortices will not be combined, RGAM = 0.0. The vortex 
reduction factor, VRF, is set to 1.0. For chine bodies the value of VRF is reset by information 

read in input Item 29(d). For chine bodies the input values of XTR1, XTR2, and EMKF are 
not used. 

The only items of interest in input Item 8 are the tolerance for the vortex trajectory 
calculation, E5, and the vortex core size, RCORE. E5 is set to 0.05, a typical value, and 
RCORE is set to 0.05. Some additional comments on the tolerance parameter E5 are 
contained in the previous section on Error Messages and Stops. 

Input Item 9 contains the flags for the source distribution calculation. NSOR = -1 
specifies the source distribution is unknown and is to be calculated by VTXCHN. NPRT = 2 
provides for additional output regarding the source distribution. This output includes a 
comparison of the input geometry and the geometry calculated from the source distribution 
which provides a good means of checking the calculated source distribution. Also printed 
are the surface velocities calculated by the source distribution. 
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Input Item 12 contains NXR, the number of entries in the body geometry table. Items 
13, 14, and 15 are the axial stations, XR, the corresponding body radii, R, and the body 
slopes, DR, respectively. The axial stations and the radii are nondimensionalized with 
respect to the body length, SL. The values of R and DR should correspond to an equivalent 
area body. 

Input Items 18(a) and 18(b) contain the mapped body radius, RO, and slope, DRO, at 
the axial stations specified in Item 13. 

Input Item 20(a) specifies that 30 mapping coefficients are used, MNFC, and the body 
cross section shape is specified at 3 axial stations, MXFC. The three axial stations where the 
cross sections are specified are read in input Item 20(b), XFC = 7.19, 13.56, 19.94. 

Input Item 21 is included for chine bodies. Three integers are specified, NCSEP, which 
indicate that the seventh corner point is taken as the chine edge for the three cross section 
shapes to be read in Items 22 and 23. 

Input Items 22 and 23 are repeated for each input cross section (MXFC = 3). Item 22, 
NR, specifies that 12 comer points are used to describe the cross section. XRC and YRC in 
Items 23(a) and 23(b) are the comer points describing the cross section. 

Input Items 28 and 29 are included for chine bodies. Item 28 contains information to 
start the chine solution. Item 28, NCSTRT, indicates that the chine initialization routine is to 
be used for the first 8 axial stations, and the boundary condition is applied at 0.25, FVCHN, 
of the exposed strake span. 

Input Item 29(a) specifies that the chine cross section starts at axial station 0.00, 
XCHNB, and that it end at axial station 30.00, XCHNE. 

Input Item 29(b) contains values of NCIR for the portion of the body in front of 
XCHNB, for the portion of the body between XCHNB and XCHNE, and for the body 


- 83 - 


sections after XCHNE. NCH (NCIR) is set to 3 for all sections of the body. The values of 
NCIR specified here overide the value specified in input Item 1. 

Input Item 29(c) specifies the boundary layer flag, NBLSEP, for the portions of the body 
before XCHNB and after XCHNE. Since the entire body has a chine cross section, these 
values are not significant. 

Input Item 29(d) specifies the vortex reduction factor, VRF, for the portions of the body 
before XCHNB and after XCHNE. Since the entire body has a chine cross section, these 
values are not significant. 

Sample Case 2 .- Chine Sample Case 2, Figure 25, is an example of a body with a 
circular nose followed by a growing chine section. 


Input Item 1 indicates that the case consists of a body with noncircular/ nonsimilar 
cross sections, NCIR = 3, and a symmetric vortex wake, ISYM = 0. The mapping flag 
NCF = 0 indicates that the numerical mapping is calculated. The separation criteria is 
laminar, NBLSEP = 1, but no reverse flow separation is to be computed, NSEPR = 0. A 
diffuse vortex core model is to be used, NDFUS = 1, and the two-dimensional unsteady 
pressure term is included in pressure calculation, NDPHI = 1. INP = 0 indicates that the 
nose forces ahead of the first axial station will be estimated. NXFV = 0 and NFV = 0 indicate 
that no special output stations or field points are to be read in Items 30 and 31. No inital 
vortices are to be read, NVP = NVR = NVM = NVA = 0, and asymmetric shedding is not 
forced, NASYM = 0. 

Input Item 2 contains two additional flags. NVPHI = -1 indicates that all vortices are 
considered independent for the d*/dt calculation, and NCHINE = 1 indicates that the body 
has a chine cross section portion. 

Input Item 3 indicates there will be five title cards input in Item 3. NPRNTP = 1 
specifies the pressure distribution be output at all x stations, and NPRNTS = 1 causes a 
summary of the separation points to be output at all x stations. The flag NPRNTV = 1 
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results in a summary of the vortex field at all axial stations. NPLOTV — 3 and NPLOTA 2 
specifies that a line printer plot of the vortex field is created at all x stations if separation has 
occurred. NTH = 0 indicates that the circle d's are calculated every 5 degrees by the 
program, and NCORE = 0 set the maximum core size to 0.05*D. Skin friction is not 
estimated, NSKIN = 0, and NCOMP = 0 specifies an incompressible calculation. 

Input Item 4 is a set of five title/summary cards which are printed at the top of the first 
output page. 

Input Item 3 contains the reference information for the run. The reference area, 
REFS = 3.14159, is the cross sectional area at the end of the circular portion of the body, the 
reference length, REFL = 2.00, is the diameter at the end of the circular portion of the body, 
and the moment center, XM, is located at station 10.00. The length of the body, SL, is 20.00. 

Item 6 contains the flow conditions. The angle of attack, ALPHAC, is 20 degrees, and 
there is no roll angle, PHI = 0.0. The Reynolds number, RE, is 1000000, and the Mach 
number is 0.0. 

Item 7 specifies run parameters. The calculation procedure begins at x station, XI, 0.50 
and continues to station, XF, 13.0 in DX increments of 0.50. A transition region is not 
specified, XTR1 = XTR2 = 0.0. Separation vortices are initially placed a minimum 1.05 radii, 
EMKF, from the axis of the missile. Vortices will not be combined, RGAM = 0.0. The vortex 
reduction factor, VRF, is set to 1.0. For the circular portion of the body the vortex reduction 
factor will be reset in input Item 29(d). 

The only items of interest in input Item 8 are the tolerance for the vortex trajectory 
calculation, E5, and the vortex core size, RCORE. E5 is set to 0.03, a typical value, and 
RCORE is set to 0.05. 

Input Item 9 contains the flags for the source distribution calculation. NSOR = -1 
specifies the source distribution is unknown and is to be calculated by VTXCHN. NPRT = 2 
provides for additional output regarding the source distribution.. 
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Input Item 12 contains NXR, the number of entries in the body geometry table. Items 
13, 14, and 15 are the axial stations, XR, the corresponding body radii, R, and the body 
slopes, DR, respectively. The axial stations and the radii are nondimensionalized with 
respect to the body length, SL. The values of R and DR should correspond to an equivalent 
area body. 

Input Items 18(a) and 18(b) contain the mapped body radius, RO, and slope, DRO, at 
the axial stations specified in Item 13. 

Input Item 20(a) specifies that 25 mapping coefficients are used, MNFC, and the body 
cross section shape is specified at 5 axial stations, MXFC. The five axial stations where the 
cross sections are specified are read in input Item 20(b), XFC = 6.00, 7.00, 13.00, 19.00, and 
20 . 00 . 

Input Item 21 is included for chine bodies. Three integers are specified, NCSEP, which 
indicate that the tenth corner point is taken as the chine edge for the five cross section 
shapes to be read in Items 22 and- 23. 

Input Items 22 and 23 are repeated for each input cross section (MXFC = 5). Item 22, 
NR, specifies that 19 comer points are used to describe the cross section. XRC and YRC in 
Items 23(a) and 23(b) are the comer points describing the cross section. 

Input Items 28 and 29 are included for chine bodies. Item 28, NCSTRT, indicates that 
the chine initialization routine is not used. 

Input Item 29(a) specifies that the chine cross section starts at axial station 6.00, 
XCHNB, and that it end at axial station 20.00, XCHNE. 

Input Item 29(b) contains values of NCIR for the portion of the body in front of 
XCHNB, for the portion of the body between XCHNB and XCHNE, and for the body 
sections after XCHNE. NCH (NCIR) is set to 0 (circular) for the portion of the body in front 
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of station XCHNB. For XCHNB < X < XCHNE NCH is set to 3, and for the portion of the 
body after XCHNE NCH is set to 0. 

Input Item 29(c) specifies the boundary layer flag, NBLSEP, for the portions of the body 
before XCHNB and after XCHNE. For both portions of the body the separation criteria is 
laminar, NCBL(l) = NCBL(3) = 1. 

Input Item 29(d) specifies the vortex reduction factor, VRF, for the portions of the body 
before XCHNB and after XCHNE. For both portions of the body the vortex reduction factor 
is VRFC(l) = VRFC(3) = 0.6. 


Output Description 

A typical chine body output file from VTXCHN is described in this section. In general, 
the output quantities from VTXCHN are labeled and each page is headed with approximate 
descriptive information. The output from Chine Sample Case 1, Figure 24, requires 151 
pages, and it will be described briefly page by page. Representative pages are shown in 
Figure 26. 

The first output page contains the title/summary cards, the reference quantities, the 
flow conditions, the initial run conditions, and the run flags. Page 1, Figure 26(a), is simply 
an echo of the input parameters. Also contained on this page are the force and moment 
definitions. Page 2, Figure 26(b), echos the input body geometry, XR, R, and DR. The 
dimensional X and R values are also printed. For noncircular bodies, R and DR, correspond 
to an equivalent area body. 

Pages 3 through 9, Figures 26(c) - 26(i), contain information from the three-dimensional 
source distribution calculation. Page 3, Figure 26(c), is a summary of the source placement 
and the body geometry at these locations. Page 4, Figure 26(d), shows the source locations 
and strengths, and pages 5 and 6, Figures 26(e) and 2h(f) are a comparison of the input 
geometry shape, R/L(INPT), and the shape calculated from the source distribution 
R/L(S.D.). This comparison is indicative of the quality of the source distribution. A second 
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indication of difficulties with the source distribution will show up in the number of 
iterations required to find the proper shape, M-INDEX. If this number is less than four or 
five, the chances are good the source distribution is a good representation of the body. If 
this number is large, examine the input and predicted shapes carefully for problems. Pages 
7 and 8, Figures 26(g) and 26(h), are the surface velocities calculated by the source 
distribution, and they should be examined to ensure that no exceedingly large velocities 
occur, another indication of a problem with the body model. Small oscillations in the 
velocities on the cylindrical portion of a body are to be expected. 

Page 9, Figure 26(i), contains a summary of source locations and strengths, and the 
values on this page are in the correct format for input into VTXCHN for additional cases 
using the same geometry and Mach number. Also shown on this page are the mapped 
radius, RO and slope DRO. 

Pages 10 through 23, Figures 26(j) - 26(w), contain output from the numerical mapping 
portion of the program. Page 10, Figure 26(j), contains the input and mapped comer points 
for the first noncircular chine cross section input, X = 7.19. Page 11, Figure 26(k), contains 
the mapping coefficients calculated, and pages 12 and 13, Figure 26(1) and 26(m), display the 
circle and body points used in the pressure calculation. For chine cross sections the chine 
edge point is indicated. Pages 14 through 21, Figures 26(n) - 26(u), contain the same 
information as pages 10 through 13, Figure 26(j) - 26(m), for the two additional chine 
sections input, X = 13.56 and 19.94. 

Pages 22 and 23, Figures 26(v) and 26(w), contain a summary of the numerical mapping 
comer points and mapping coefficients. The information on these two pages can be used as 
input to the code if the same body needs to be run for different flow conditions. 

Page 24, Figure 26(x), contains a summary of chine information. Information on this 
page describes how the chine solution is started, beginning and end of chine section, and 
program control variables for the different parts of the body. 
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Page 25, Figure 26(y), shows the pressure distribution at the first axial station before a 
vortex has been placed to satisfy the Kutta condition at the chine edge. The vortex required 
to satisfy the Kutta condition is printed at the bottom of this page. 

Page 26, Figure 26(z), starts the output from the analysis part of the code. The first 
numeric line on page 26 contains the axial locations, the local radius, the body slope, and the 
transition Reynolds number, RETR. The predicted results on page 26 are a summary of the 
velocities and pressures on the circumference at this particular x station. The coordinates 
are dimensional, and the velocities are made dimensionless by the freestream velocity. The 
pressure coefficient is followed by the unsteady term d$/dt, velocities in the circle plane, 
and the value of the numerical mapping CS. The sharp chine edge is indicated. 

Following the pressure distribution is a summary of the section loadings and the 
overall loadings up to the current section. CN(X), CY(X), and CA(X) are the normah, side-, 
and axial-force coefficients per unit length, respectively. CN, CY, and CA are the overall 
normal-, side-, and axial-force coefficients from the nose to the current x station, 
respectively. CM, CR, and CSL are the pitching-, yawing-, and rolling-moment coefficients, 
respectively. XCPN and XCPY are the x locations of the centers of pressures for the normal 
and side forces, respectively. 

Page 27, Figure 26(aa), contains a line printer plot of the vortex field at the first x 
station. The tracking of the vortex to the next axial station takes place after this line printer 
plot. 


Page 28, Figure 26(ab), a summary of the vortex field at the next axial station, X = 1.0. 
The new vortex required to satisfy the Kutta condition at this station follows the summary 
of the vortex field table. 

Pages 28-117 contain output for axial stations 1.0 to 14.5. The output for each axial 
station has the same form as that described for axial station 0.5 on pages 26 through 28, 
Figures 26(z) - 26(ab). The output for axial stations 4.5 and 14.5 are shown on pages 50 - 53, 
Figure 26(ac) - 26(af), and 114-117, Figures 26(ag) - 26(aj), respectively. 
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Page 118, Figure 26(ak), of the output file depicts the vortex field at the final axial 
station and contains a summary of the overall loads. Page 119, Figure 26(al), shows the 
pressure distribution at the final output station. This pressure distribution is required if a 
restart calculation is to be run starting at this final station. 

CONCLUSIONS 

An engineering prediction method based on a rational flow modeling technique and 
the associated code VTXCHN to predict the vortex shedding from forebodies with chines in 
subsonic flow at angles of attack and roll are described. Comparisons of measured and 
predicted aerodynamic characteristics and flow field quantities are used to verify the flow 
model and prediction method for a variety of configurations under a wide range of flow 
conditions. The method has proved successful in representing the principal features of the 
complex flow field on the lee side of these configurations at moderate angles of attack; 
therefore, it has application as an engineering or preliminary design technique directed at 
the prediction of nonlinear aerodynamic characteristics resulting from high angle of attack 
flows. 

The prediction method described herein has further application as one component of a 
larger prediction method for complete configurations consisting of an arbitrary body and 
multiple fins. The ability to model the correct flow field in the vicinity of the body leads to 
the capability to calculate body separation wake-induced interference effects on fins and 
other control surfaces. The vortex shedding analysis and prediction techniques are also 
applicable for use in other codes; for example, the methods developed in this investigation 
can be transferred to higher order codes such as panel codes. This was demonstrated in the 
vortex shedding work reported in Reference 3 and 16. 

Finally, the prediction method presented in this report must be considered as 
preliminary. Even though verification of the method by comparisons with experiment has 
demonstrated some success, there are configurations and flow conditions which must be 
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approached with caution. For this reason, several recommendations for improvements to 
the method are suggested in the following section. 


RECOMMENDATIONS 

In the course of development and verification of the code VTXCHN, several specific 
areas requiring additional work were identified. Since the additional effort was beyond the 
scope of the present investigation, they are noted below as recommendations for future 
work. 

The first and most important recommendation is for a thorough testing and further 
verification of VTXCHN. The code should be applied to a variety of configurations and 
flow conditions for which experimental results are available for comparison purposes. This 
should better define the operational limits of the method. The major difficulty will be the 
lack of useful experimental results on appropriate forebodv configurations. 

The numerical difficulty associated with secondary separation discussed in a previous 
section needs further investigation. A minor modification may correct the problem, but 
some additional effort is required to understand the source of the problem. 

The prediction method may have additional uses which have not been explored. For 
example, the code has vortex shedding technology built into it which can be transferred to 
three-dimensional panel codes. The code also has potential use in providing separation line 
and vortex wake information for starting solutions for higher order prediction methods. 

Finally, the prediction method may have application in the investigation of other 
complex vortex phenomena such as asymmetry and vortex bursting. 
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Figure 1 Lee side vortex formation on a forebody with chine. 
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(b) Numerical transformation procedure 
Figure 2.- Conformal mapping nomenclature. 
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Figure 7.- Vortex cloud modification for d</>/dt calculation. 
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Figure 1 9.- Subroutine calling sequence for program VTXCHN. 
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Figure 20.- Subroutine cross reference list for program VTXCHN. 
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Figure 20.- Continued 
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Figure 20.- Concluded. 
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Figure 21.- Common block cross reference list for program 
VTXCHN. 
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Figure 21.- Continued. 
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Figure 21.- Concluded. 
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Figure 22.- Input forms for program VTXCHN. 
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Figure 22.- Continued. 
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Figure 22.- Concluded. 
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0.00586 

0.23980 

0.00000 

-0.07482 

-0.28713 


0.0 

0.0 


0.07834 

0.34060 

1.08692 

1.27166 

0.02061 

0.03799 

0.03463 

. WvW 

0.20492 

0.00000 

-0.06198 

-0.31992 

0.02069 

0.04134 

0.03715 

V • WWV 
0.21662 
/a 6UJUH0UH 

V • VVX/W 

-0.09862 

-0.31156 


1.06 

0.0 


0.09796 

0.47684 

1.13379 

0.02419 

0.03799 

0.03069 

0.18673 

« AOLOUXOk 
V ■ X/xPOvv 

-0.11766 

0.02461 

0.04134 

0.03146 

0.19376 

0.00000 

-0.14810 


0.0 

0.0 


0.11744 

0.61308 

1.15422 

0.02791 

0.03799 

0.02791 

0.16184 

0.00000 

-0.16184 

0.02817 

0.04134 

0.02817 

0.17136 

0.00000 

-0.17135 


1.00 

0.05 


0.13787 

0.74932 

1.17371 

0.03059 

0.03799 

0.02419 

0.11766 

0.00000 

-0.18673 

0.03146 

0.04134 

0.02460 

0.14810 

0.00000 

-0.19377 


64333E+00 1.96750E+00 
29000E+00 0 . 00000 E-c 00 
76125E+00-0 . 29725E+00 
72792E+00 1 . 83667E+00 

35626E+00 2 . 90000E+00 
39875E+00 0.00000E+00 
13583E+00-0 . 43500E+00 
51333E+00 2.69792E+00 

41107E+00 3 . 16905E+00 
46060E+00 0 .00000E+00 
23073E+00-0.48274E+00 
, 61210E+00 2 . 65240E+00 


0.00 


30.00 

3 

3 

3 

1 

1.00 

1 

1.00 


Figure 24.- Sample case 1 input file. 
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INPUT ITEM 


(1) 300101 10 

(2) -110 

(3) 3 1 1 1 3 2 0 0 

(4) NIELSEN ENGINEERING k RESEARCH, INC. 

CIRCULAR - CHINE BODY 
SAMPLE CASE 2 

(6) 3.14169 

(6) 20.000 

(7) 0.60 

(8) .030 

(9) -1 2 

(12) 36 

(13) 0.00000 
0.09926 
0.35 
1.16000 
1.29344 

(14) 0.00000 
0.02796 

.05638 
0.03776 
0.00222 

(16) 0.34286 

0.22231 
.02600 
-0.16434 
-0.33452 

(18a) 0.00000 

0.02796 
.05209 
0.03776 
0.00222 

(18b) 0.34286 

0.22231 
.01640 
-0.16434 
-0.33452 
26 5 

6.00 
10 10 
19 

S*??? 00 !* 00 0 . 17365E+00 0 . 34202E+00 0.60000E+00 0.64279E+00 0 
0.86603E+00 0.93969E+00 0.98481E+00 0.10000E+01 0.98481E+00 0 

S'SSS 03 !* 00 0 - 7 ®®04E^00 0.64279E+00 0.60000E+00 0.34202E+00 0 
0 . 00000E+00 

(23b) -0 . 10000E+01-0 . 98481E+00-0 . 93969E+00-0 . 86603E+00-0 . 76604E+00-0 
-0.60000E+00-0.34202E+00-0.17366E+00 0.00000E+00 0.17366E+00 0 
0 . 60000E+00 0.64279E+00 0.76604E+00 0.86603E+00 0.93969E+00 0. 
0.10000E+01 

(22) 19 

(23a) 0 . 00000E+00 0.17365E+00 0 . 34202E+00 0.50000E+00 0.64279E+00 0 

0 . 866C3E+00 0 . 93969E+00 0.98481E+00 0.13000E+01 0.98481E+00 0. 
0 . 86603E+00 0 . 76604E+00 0.64279E+00 0 . 50000E+00 0.34202E+00 0 
0 . 00000E+00 

(23b) -0 . 10000E+01-0 . 98481E+00-0 . 93969E+00-0 . 86603E+00-0 . 76604E+00-0 . 
-0 . 50000E+00-0 . 34202E+00-0 . 17365E+00 0.00000E+00 0. 17365E+00 0. 
0 . 50000E+00 0 . 64279E+00 0.76604E+00 0 . 86603E+00 0.93969E+00 0 
0 . 10000E+01 


0 0 0 
PROGRAM VTXCHN 


2.0000 

10.000 

20.000 

O 

£ . WW 




0.000 

1000000 . 

0.00 





13.000 

0.50 

0.0 

0.0 

1.05 

0.0 

1.00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.05 

0.00164 

0.00656 

0.01468 

0.02594 

0.04019 

0.06729 

0.07706 

0.12366 

0.16000 

0.17798 

0.20729 

0.23763 

0.26864 

0 . 30000 

0.66 

0.95 

1 . WVVv 

1.03136 

1.06237 

1.09271 

1.12202 

1.17634 

1 . 20074 

1.22294 

1.24271 

1.26981 

1.27406 

1.28632 

1.29836 

1.30000 






0.00056 

0.00222 

0.00490 

0.00847 

0.01276 

0.01769 

0.02273 

0.03304 

0.03776 

0.04192 

0.04634 

0.04789 

0.04947 

0 . 06000 

.05907 

.06638 

0.06000 

0.04947 

0.04789 

0.04634 

0.04192 

0.03304 

0.02796 

0.02273 

0.01769 

0.01276 

0.00847 

0.00490 

0.00056 

0.00000 






0.34076 

0.33452 

0.32426 

0.31021 

0.29265 

0.27191 

0.24834 

0.19419 

0.16434 

0.13308 

0.10073 

0.06758 

0.03392 

0 . 00000 

.00000 

-.02600 

V • WvVv 

-0.03392 

-0.06768 

-0.10073 

-0.13308 

-0.19419 

-0.22231 

-0.24834 

-0.27191 

-0.29266 

-0.31021 

-0.32426 

-0.34076 

-0.34286 






0.00066 

0.00222 

0.00490 

0.00847 

0.01276 

0.01769 

0.02273 

0.03304 

0.03776 

0.04192 

0.04634 

0.04789 

0.04947 

0 . 05000 

.06618 

.06209 

0 . 06000 

0.04947 

0.04789 

0.04634 

0.04192 

0.03304 

0.02796 

0.02273 

0.01769 

0.01276 

0.00847 

0.00490 

0.00056 

0.00000 






0.34076 

0.33452 

0.32426 

0.31021 

0.29266 

0.27191 

0.24834 

0.19419 

0.16434 

0.13308 

0.10073 

0.06758 

0.03392 

0 . 00000 

.00000 

-.01640 

0.00000 

-0.03392 

-0.06758 

-0.10073 

-0.13308 

-0.19419 

-0.22231 

-0.24834 

-0.27191 

-0.29266 

-0.31021 

-0.32426 

-0.34076 

-0.34286 






7.00 

13.00 

19.00 

20.00 




10 10 

10 







76604E+00 

93969E+00 

17366E+00 

64279E+00 

34202E+00 

98481E+00 


76604E+00 

93969E+00 

17366E+00 

64279E+00 

34202E+00 

98481E+00 


(aj Pane 1 


Figure 25.- Sample case 2 input file. 
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( 22 ) 

(23.) 


(23b) - 


( 22 ) 

(23.) 


(23b) 


( 22 ) 

(23.) 


(23b) 


(28) 

(29.) 

(29b) 

(29c) 

(29d) 


Z 00000E+00 0 . 17366E+00 0.34202E.00 0.60000E+00 0.84279E+00 0 . 76604E+00 
088803E+00 0.93969E+00 0.98481E+00 0. 18000E-.01 0.98481E+00 0.93989E-.00 
0!86603E+00 0. 76804E+00 0.84279E+00 0.60000E+00 0.34202E+00 0.17366E+00 
0 ! 00000E+00 

0 10000E-.01-0 . 98481E+00-0 . 93969E+00-0 . 88803E+00-0 . 76604E*00-0 . 64279E+00 
0 ' 60000E+00-0. 34202E+00-0. 17366E+00 0.00000E+00 0.1736SE+00 0.34202E+00 
0!b0000E+00 0.64279E+00 0.78804E+00 0.88003E-.00 0.93969E+00 0.98481E-.00 
0 . 10000E+01 
19 

0 00000E+00 0. 17366E+00 0 . 34202E+00 0.60000E+00 0.64279E+00 0.76604E+00 
0 86603E-.00 0.93969E+00 0.98481E-.00 0.13000E-.01 0.98481E+00 0.93969E+00 
0 ! 86603E+00 0.76804E+00 0.64279E+00 0.60000E+00 0.34202E+00 0.17366E+00 
0.00000E+00 

0 10000E+01-0 . 98481E+00-0 . 93969E+00-0 . 86603E+00-0 . 78604 E+00-0 .642 79E+00 
0 60000E+00-0 • 34202E-.00-0 . 1736BE+00 0.00000E+00 0.17366E+00 0.34202E+00 
0 . 60000E+00 0.64279E*00 0.78804E+00 0.86603E+00 0.93969E+00 0.98481E-.00 
0 . 10000E+01 
19 

0 00000E+00 0.17366E+00 0. 34202E+00 0.60000E+00 0.64279E+00 0.76604E+00 
0 86603E+00 0 . 93969E-C00 0.98481E+00 0 . 10000E+01 0.98481E+00 0.93969E+00 
0 . 86603E+00 0.76604E+00 0.64279E+00 0 . 60000E+00 0.34202E+00 0.17366E+00 
0 . 00000E+00 

0 . 10000E+01-0 . 98481E+00-0 . 93969E-.00-0 . 86803E+00-0. 76604E+00-0 . 64279E+00 
•0 B0000E+00-0.34202E+00-0.1738BE+00 0.00000E-.00 0. 17386E+00 0.34202E*00 
0 . 60000E+00 0 . 64279E.00 0.76604E+00 0.86603E+00 0.93969E*00 0.98481E+00 
0 . 10000E+01 
0 0 . 6000 


6.00 


20.00 

0 

3 

0 

1 

0.60 

1 

0.60 


(b) Page 2 

Figure 25 - Concluded. 
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PROGRAM VTXCHN 


NIELSEN ENGINEERING A RESEARCH, INC. 

CHINE BODY 

REF. AIAA 87-2617 

CHINE SAMPLE CASE 1 


REFERENCE QUANTITIES 


REF. AREA 22.16 

REF. LENGTH 8.70 

X (MOMENT CENTER) . 23.00 

BODY LENGTH (L) . . . 73.40 

BASE DIAMETER (D) . 8.70 

STRATFORD (LAM)... 0.087 

STRATFORD (TURB) . . 0.360 


FLOW CONDITIONS 

ALPHA C (DEG.) ... 20.00 

PHI (DEG.) 0.00 

MACH NUMBER 0.00 

REYNOLDS NO. (D) . . 1.000E+06 

VORTEX FACTOR 1.00 

ALPHA (DEG.) 20.00 

BETA (DEG.) 0.00 

U/V0 0.940 

V0*DT 0.632 

RE/D 1 . 149E+05 

RE/D SIN (ALPHA) .. 3.931E+04 


INITIAL 

CONDITIONS 

XI 

0.600 

XF 

14.600 

DX 

0.600 

XTRl 

16.960 

XTR2 

24.660 

TOL 

0.05000 

OPTIONS . 

. . 








NCIR NCF 

3 0 

ISYM 

0 

NBLSEP 

1 

NSEPR 

0 

NDFUS 

1 

NDPHI INP 

1 0 

NXFV 

0 

NFV NVP 

0 0 

NVR 

0 

NPRNTP NPRNTS 
1 1 

NPRNTV 

1 

NPLOTV 

3 

NPLOTA 

2 

NTH 

0 

NCORE NSKIN 
0 0 

NCOMP 

0 

NDADT 

0 



EMKF RGAM 

1.050 0.000 


NVM NVA NASYM 
0 0 0 


DEFINITIONS 

OF OUTPUT QUANTITIES ... 


CA (X) 

(DCA/DX) / (Q.O) , POSITIVE AFT 


CN (X) 

(DCN/DX) / (Q.D; , POSITIVE UP 


CY (X) 

(DCY/DX)/(Q»D) , POSITIVE RIGHT 


CA 

FX/ (Q.REFS) , POSITIVE AFT 


CN 

FZ/ (Q.REFS) , POSITIVE UP 


CY 

FY/ (Q.REFS) , POSITIVE RIGHT 


CM 

M/ (Q.REFS.REFL) , POSITIVE NOSE 

UP 

CR 

R/ (Q.REFS.REFL) , POSITIVE NOSE 

RIGHT 

CLL 

L/ (Q.REFS.REFL) , POSITIVE COUNTERCLOCKWISE 


(a) Page 1 


RCORE/D 

0.050 


NVPHI NCHINE 
-1 1 


Figure 26.- Sample case 1 output. 
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INPUT GEOMETRY CHARACTERISTICS 
NXR = 29 


X/L 

REQ/L 

DREQ/DX 

X 

REQ 

0 . 0000 

0.0000 

0.3199 

0.0000 

0.0000 

0.0196 

0 . 0060 

0.2944 

1.4401 

0.4426 

0.0391 

0.0115 

0.2591 

2.8699 

0.8448 

0.0687 

0.0161 

0.2318 

4.3100 

1.1864 

0.0783 

0.0206 

0.2049 

6.7502 

1.6128 

0.0980 

0.0242 

0.1867 

7.1903 

1.7755 

0.1174 

0.0279 

0.1618 

8.6201 

2.0486 

0.1379 

0.0306 

0.1176 

10.1197 

2.2463 

0.1847 

0.0346 

0.0620 

13.5599 

2.6418 

0.2035 

0.0356 

0.0471 

14.9398 

2.6138 

0.2350 

0.0368 

0.0354 

17.2497 

2.7026 

0.2717 

0.0380 

0.0207 

19.9398 

2.7885 

0.3406 

0.0380 

0.0000 

25.0000 

2.7886 

0.4768 

0.0380 

0.0000 

35.0001 

2.7886 

0.6131 

0.0380 

0.0000 

46.0001 

2.7886 

0.7493 

0.0380 

0.0000 

66.0001 

2.7886 

0.9312 

0.0380 

0.0000 

68.3601 

2.7885 

1 . 0000 

0.0380 

-0.0207 

73.4000 

2.7885 

1.0366 

0.0368 

-0.0354 

76.0901 

2.7026 

1.0681 

0.0356 

-0.0471 

78.4000 

2.6138 

1.0869 

0.0346 

-0.0620 

79.7799 

2.5418 

1.1338 

0.0306 

-0.1176 

83.2202 

2.2463 

1.1542 

0.0279 

-0.1618 

84.7197 

2.0486 

1.1737 

0.0242 

-0.1867 

86.1603 

1.7766 

1.1933 

0.0206 

-0.2049 

87.6897 

1.6128 

1.2129 

0.0161 

-0.2318 

89.0298 

1.1864 

1.2326 

0.0115 

-0.2591 

90.4699 

0.8448 

1.2520 

0 . 0060 

-0.2944 

91.8997 

0.4426 

1.2717 

0 . 0000 

-0.3199 

93.3398 

V . WVKf 


(b) Page 2 


Figure 26.- Continued. 
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SO'-RCE LOCATIONS AND BODY RADIUS AND SURFACE SLOPE AT THESE LOCATIONS 
SOURCE SPACING IS 0.76000 TIMES LOCAL RADIUS 


X/L 

0.00200 

0.00246 

0.00303 

0.00373 

0.00459 

0.00664 

0.00694 

R/L 

0.00061 

0.00076 

0.00093 

0.00115 

0.00141 

0.00173 

0.00213 

DR/DX 

0.31732 

0.31672 

0.31698 

0.31507 

0.31396 

0.31268 

0.31089 

X/L 

0.00864 

0.01061 

0.01294 

0.01692 

0.01969 

0.02410 

0.02967 

R/L 

0.00263 

0.00323 

0.00398 

0.00489 

0.00602 

0.00729 

0.00883 

DR/DX 

0.30881 

0.30625 

0.30309 

0.29922 

0.29444 

0.28627 

0.27636 

X/L 

0.03619 

0.04421 

0.06375 

0.06498 

0.07816 

0.09368 

0.11113 

R/L 

0.01069 

0.01272 

0.01497 

0.01757 

0.02067 

0.02339 

0.02670 

DR/DX 

0.26434 

0.25196 

0.23872 

0.22323 

0.20617 

0.19079 

0.16991 

X/L 

0.13116 

0.16344 

0.17739 

0.20288 

0.22956 

0.25702 

0.28516 

R/L 

0.02971 

0.03193 

0.03400 

0.03668 

0.03661 

0.03752 

0.03799 

DR/DX 

0.13211 

0.09910 

0.07070 

0.04764 

0.03740 

0.02668 

0.01668 

X/L 

0.31366 

0.34215 

0.37064 

0.39913 

0.42763 

0.46612 

0.48461 

R/L 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

DR/DX 

0.00811 

0.00000 

0.00000 

0.00000 

0.00000 

Ok OUAOkOkOk 

0.00000 

X/L 

0.61310 

0.64160 

0.67009 

0.69858 

0.62707 

0.65557 

0.68406 

R/L 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

DR/DX 

0.00000 

0.00000 

0.00000 

0.00000 

V . KWvfcW 

0.00000 

0.00000 

X/L 

0.71266 

0.74104 

0.76954 

0.79803 

0.82652 

0.86601 

0.88351 

R/L 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

DR/DX 

0.00000 

ot 

V • PlfWvtf 

0.00000 

V . Wvvv 

0.00000 

0.00000 

0.00000 

X/L 

0.91200 

0.94049 

0.96898 

0.99748 

1.02697 

1.05384 

1.08096 

R/L 

0.03799 

0.03799 

0.03799 

0.03799 

0.03716 

0.03616 

0.03494 

DR/DX 

0.00000 

-0.00280 

-0.01138 

-0.01996 

-0.03110 

-0.04179 

-0.06727 

X/L 

1.10717 

1.13183 

1.16490 

1.17673 

1.19360 

1.20901 

1.22179 

R/L 

0.03288 

0.03076 

0.02778 

0.02382 

0.02056 

0.01704 

0.01406 

DR/DX 

-0.08699 

-0.11623 

-0.16271 

-0.18861 

-0.20630 

-0.22642 

-0.24410 

X/L 

1.23233 

1.24101 

1.24786 

1.26326 

1.25750 

1.26077 

1.26328 

R/L 

0.01166 

0.00913 

0.00721 

0.00566 

0.00435 

0.00336 

0.00258 

DR/DX 

-0.26876 

-0.27438 

-0.28681 

-0.29699 

-0.30151 

-0.30676 

-0.30902 

X/L 

1.26621 







R/L 

0.00198 







DR/DX 

-0.31163 








X (ST) 

XS(1) 

X (RM) 

R (MAX) 

XS(L) 

X (ST) 



0.00169 0.00200 0 

.27166 0 

.03799 1 

.26666 1 

.26620 



(c) Page 3 


Figure 26.- Continued. 
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FOR THIS CASE THERE ARE 71 SOURCES 


X/L 

Q 

2 . 0000E-03 
6 . 9391E-08 

2 . 4610E-03 
1 . 5171E-08 

3 . 0283E-03 
8 . 7780E-08 

3 . 7263E-03 
9 . 3363E-08 

X/L 

Q 

6 . 9427E-03 
3 . 7472E-07 

8.6430E-03 
B .6B27E-07 

1.0512E-02 

8.6412E-07 

1 . 293EE-02 
1.31E2E-08 

X/L 

Q 

2 . 4101E-02 
4 . 6911E-08 

2 . 9569E-02 
6.B068E-08 

3.6190E-02 

9.4676E-00 

4.4209E-02 

1.3641E-06 

X/L 

Q 

7 . 81S7E-02 
3 . 0421E-06 

9.3S83E-02 
3 . 8366E-06 

1 . 1113E-01 
4 . 8201E-0S 

1 . 3116E-01 
E . 9922E-0E 

X/L 

Q 

2 . 0288E-01 
3 . 4914E-08 

2.29B6E-01 

2.7479E-06 

2 . 6702E-01 
1 . 3126E-06 

2.8E16E-01 

6.36B6E-06 

X/L 

Q 

3 . 7064E-01 
4.1969E-06 

3 . 9913E-01 
-6.1886E-06 

4 . 2763E-01 4 . B612E-01 

4 . 3806E-00 -4 . 9101E-00 

X/L 

Q 

E.4160E-01 
4 . 6168E-00 

6 . 7009E-01 
“4 . 7422E-06 

B . 98S8E-01 6 . 2707E-01 

4 . 8063E-00 -4.6986E-06 

X/L 

Q 

7 . 12B5E-01 
4 . 761SE-06 

7.4104E-01 

-4.6042E-06 

7 . 0954E-01 7 . 9803E-01 

4.8172E-06 -4 . B173E-06 

X/L 

Q 

8 . 8361E-01 
6 . 2932E-00 

9 . 1200E-01 
-4.3246E-06 

9 . 4049E-01 9 . 6898E-01 

1 . 2747E-0B -1 . 8080E-05 

X/L 

« 

1 . 0538E+00 
-9 . 2092E-06 

1.0810E+00 
-2 . 7861E-06 

1 . 1072E+00 1 . 1318E+00 

-3.0873E-06 -4 .6932E-06 

X/L 

Q 

1 .1936E+00 
-3 . 6964E-06 

1.2090E+00 

-2.6898E-06 

1 . 2218E+00 1 . 2323E+00 

-1 .8986E-06 -1 . 3670E-06 

X/L 

1 .2633E+00 

1 . 257SE+00 

1 . 2608E+00 

1.2633E+00 


Q -3.438SE-06 -1.90S6E-06 -1.3171E-0S -4.S432E-07 


4.E8E2E-03 6.6422E-03 

1.6719E-07 2.3969E-07 

1.S917E-02 1 . 9686E-02 

1 . 9884E-06 3 . 1042E-06 

S.3747E-02 8.4978E-02 

1 . 7S09E-0S 2.4069E-06 

1 . S344E-01 1 . 7739E-01 

2.0738E-06 5.0476E-0S 

3 . 1386E-01 3.421SE-01 

7.3676E-08 -9.1E44E-06 

4 . 8461E-01 6 . 1310E-01 

4 . 6378E-08 -4.8024E-06 

6 . 6E67E-01 6.8406E-01 

4 . 7063E-00 -4.6E70E-06 

8 . 26S2E-01 8.ES01E-01 

4 . 9399E-00 -4 . 33E1E-08 

9 . 9748E-01 1 . 0260E+00 

2 . 8700E-08 -2.9473E-0E 

1 . 1649E+00 1 . 17E7E+00 

-8 . 2039E-0E -4.4803E-0S 

1 . 24 10E+00 1 . 2479E+00 

-8 . 60E8E-08 -E . EE27E-06 

1 . 20E2E+00 
-7 .6060E-07 
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SHAPE CALCULATED FROM SOURCE DISTRIBUTION 


X/L 

0.01000 

0.02000 

0 . 03000 

0.04000 

0.06000 

0.06000 

0.07000 

R/L(S.D.) 

0.00307 

0.00608 

0.00894 

0.01162 

0.01413 

0.01649 

0.01872 

R/L (INPT) 

0.00307 

0.00614 

0.00895 

0.01172 

0.01409 

0.01644 

0.01871 

M- INDEX 

1 

2 

1 

2 

2 

2 

1 

X/L 

0.08000 

0.09000 

0.10000 

0.11000 

0.12000 

0.13000 

0.14000 

R/L (S .D . ) 

0.02082 

0.02281 

0.02470 

0.02647 

0.02811 

0.02969 

0.03086 

R/L (INPT) 

0.02091 

0.02274 

0.02468 

0.02649 

0.02826 

0.02966 

0.03077 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.15000 

0.16000 

0.17000 

0.18000 

0.19000 

0.20000 

0.21000 

R/L (S . D . ) 

0.03193 

0.03286 

0.03370 

0.03446 

0.03610 

0.03662 

0.03606 

R/L (INPT) 

0.03164 

0.03260 

0.03336 

0.03422 

0.03490 

0.03643 

0.03686 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.22000 

0 . 23000 

0 . 24000 

0.26000 

0.26000 

0.27000 

0 . 28000 

R/L (S . D . ) 

0.03645 

0 . 03684 

0.03720 

0.03762 

0.03780 

0.03804 

0.03823 

R/L (INPT) 

0.03624 

0.03663 

0.03698 

0.03730 

0.03762 

0.03794 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.29000 

0 . 30000 

0.31000 

0.32000 

0.33000 

0.34000 

0.36000 

R/L (S .D . ) 

0.03839 

0.03853 

0.03864 

0.03872 

0.03876 

0.03878 

0.03879 

R/L (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.36000 

0.37000 

0.38000 

0.39000 

0.40000 

0.41000 

0.42000 

R/L (S . D . ) 

0.03879 

0.03879 

0.03880 

0.03880 

0.03879 

0.03879 

0.03879 

R/L (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.43000 

0.44000 

0.46000 

0.46000 

0.47000 

0.48000 

0.49000 

R/L (S . D . ) 

0.03880 

0.03880 

0.03880 

0.03879 

0.03879 

0.03879 

0.03880 

R/L (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.60000 

0.61000 

0.62000 

0.63000 

0.64000 

0.66000 

0.66000 

R/L($.D.) 

0.03880 

0.03880 

0.03879 

0.03879 

0.03880 

0.03880 

0.03880 

R/L (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.57000 

0.58000 

0.69000 

0.60000 

0.61000 

0.62000 

0.63000 

R/L (S .D . ) 

0.03880 

0.03879 

0.03879 

0.03880 

0.03880 

0.03880 

0.03879 

R/L (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.64000 

0.65000 

0.66000 

0.67000 

0.68000 

0.69000 

0.70000 

R/L (S . D . ) 

0.03879 

0.03879 

0.03880 

0.03880 

0.03880 

0.03879 

0.03879 

R/L (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0.71000 

0 . 72000 

0.73000 

0.74000 

0.75000 

0.76000 

0.77000 

R/L (S . D . ) 

0.03879 

0.03880 

0.03880 

0.03880 

0.03879 

0.03879 

0.03880 

R/L (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 

X/L 

0 . 78000 

0.79000 

0 . 80000 

0.81000 

0.82000 

0.83000 

0.84000 

R/L (S . D . ) 

0.03880 

0.03880 

0.03879 

0.03879 

0.03879 

0.03880 

0.03880 

R/L. (INPT) 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

0.03799 

M-INDEX 

2 

2 

2 

2 

2 

2 

2 
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X/L 

R/US.D.) 

R/L (INPT) 
M- INDEX 

0.85000 

0.03880 

0.03799 

2 

0.86000 

0.03879 

0.03799 

2 

0.87000 

0.03879 

0.03799 

2 

0 . 88000 
0.03879 
0.03799 
2 

0.89000 

0.03880 

0.03799 

2 

0.90000 

0.03880 

0.03799 

2 

0.91000 

0.03879 

0.03799 

2 

X/ L 

R/L (S .D . ) 
R/L (INPT) 
M- INDEX 

0.92000 

0.03879 

0.03799 

2 

0.93000 

0.03879 

0.03799 

2 

0.94000 

0.03878 

0.03799 

2 

0.95000 

0.03876 

0.03799 

2 

0.96000 

0.03868 

0.03799 

2 

0.97000 

0.03866 

0.03799 

2 

0.98000 

0.03842 

0.03799 

2 

X/L 

R/L (S .D . ) 
R/L (INPT) 
M-INDEX 

0.99000 

0.03828 

0.03799 

2 

1.00000 

0.03808 

0.03799 

2 

1.01000 

0.03787 

0.03767 

2 

1.02000 

0.03760 

0.03735 

2 

1.03000 

0.03727 

0.03703 

2 

1.04000 

0.03692 

0.03669 

2 

1.06000 

0.03664 

0.03631 

2 

X/L 

R/L (S .D . ) 
R/L (INPT) 
M-INDEX 

1 . 06000 
0.03614 
0.03592 
2 

1.07000 

0.03569 

0.03551 

2 

1.08000 

0.03617 

0.03499 

2 

1.09000 

0.03465 

0.03436 

2 

1 . 10000 
0.03384 
0.03360 
2 

1.11000 

0.03302 

0.03264 

2 

1 . 12000 
0.03208 
0.03178 
2 

X/L 

R/L(S.D ' 
R/L (INPT) 
M-INDEX 

1.13000 

0.03101 

0.03092 

2 

1.14000 

0.02977 

0.02978 

1 

1.15000 
0 . 02836 
0.02846 
2 

1 . 16000 
0 . 02674 
0.02681 
2 

1.17000 

0.02498 

0.02490 

2 

1.18000 

0.02311 

0.02304 

2 

1.19000 

0.02114 

0.02122 

2 

X/L 

R/L (S .D . ) 
R/L (INPT) 
M-INDEX 

1.20000 

0.01905 

0.01909 

2 

1.21000 

0.01684 

0.01682 

1 

1 . 22000 
0.01449 
0.01448 

2 

1.23000 

0.01199 

0.01212 

2 

1.24000 

0.00932 

0.00942 

2 

1.26000 

0.00646 

0.00660 

2 

1.26000 

0.00335 

0.00358 

2 
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VELOCITIES ON THE BODY SURFACE 

VELOCITIES 




AXIAL 

RADIAL 

TANGEN 

NORMAL 

X/L 

R/L 

VA/V0 

VR/V0 

VT/V0 

VN/V0 

0.00169 

0.00062 

-0.33774 

0.20448 

-0.25997 

0.29715 

0.01000 

0.00307 

-0.18808 

0.24935 

-0.10664 

0.29366 

0.02000 

0.00614 

-0.14678 

0.24948 

-0.07052 

0.28073 

0.03000 

0.00895 

-0.11708 

0.24305 

-0.04830 

0.26643 

0.04000 

0.01172 

-0.09406 

0.23346 

-0.03280 

0.24954 

0 . 06000 

0.01409 

-0.07646 

0.22568 

-0.01984 

0.23704 

0 . 06000 

0.01644 

-0.06090 

0.21620 

-0.01088 

0.22435 

0.07000 

0.01871 

-0.04709 

0.20610 

-0.00244 

0.21140 

0.08000 

0.02091 

-0.03612 

0.19605 

0.00465 

0.19912 

0 . 09000 

0.02274 

-0.02496 

0.18963 

0.01163 

0.19091 

0.10000 

0.02458 

-0.01360 

0.18132 

0.01946 

0.18078 

0.11000 

0.02649 

-0.00109 

0.17101 

0.02781 

0.16873 

0.12000 

0.02825 

0.01174 

0.16729 

0.03689 

0.16369 

0.13000 

0.02956 

0.02429 

0.14177 

0.04298 

0.13726 

0.14000 

0.03077 

0.03174 

0.12047 

0.04630 

0.11606 

0.16000 

0.03164 

0.03170 

0.10280 

0.04208 

0.09900 

0.16000 

0.03250 

0.02952 

0.09166 

0.03772 

0.08849 

0.17000 

0.03336 

0.03040 

0.08398 

0.03696 

0.08131 

0.18000 

0.03422 

0.03416 

0.07363 

0.03905 

0.07116 

0.19000 

0.03490 

0.03603 

0.06992 

0.03843 

0.06780 

0 . 20000 

0.03643 

0.03104 

0.04834 

0.03341 

0.04673 

0.21000 

0.03586 

0.02648 

0.04219 

0.02734 

0.04101 

0.22000 

0.03624 

0.02206 

0.04026 

0.02369 

0.03931 

0.23000 

0.03663 

0.02186 

0.03876 

0.02328 

0.03791 

0.24000 

0.03698 

0.02279 

0.03640 

0.02396 

0.03462 

0.25000 

0.03730 

0.02313 

0.03096 

0.02403 

0.03027 

0.26000 

0.03762 

0.02286 

0.02634 

0.02362 

0.02676 

0.27000 

0.03794 

0.02198 

0.02188 

0.02244 

0.02140 

0.28000 

0.03799 

0.02074 

0.01815 

0.02107 

0.01777 

0.29000 

0.03799 

0.01952 

0.01510 

0.01976 

0.01481 

0 . 30000 

0.03799 

0.01856 

0.01247 

0.01870 

0.01224 

0.31000 

0.03799 

0.01795 

0.00983 

0.01804 

0.00966 

0.32000 

0.03799 

0.01733 

0.00666 

0.01737 

0.00664 

0.33000 

0.03799 

0.01686 

0.00318 

0.01587 

0.00313 

0 . 34000 

0.03799 

0.01329 

0.00060 

0.01329 

0 . 191719010 

0.36000 

0.03799 

0.01064 

-0.00022 

0.01054 

-0.00022 

0.36000 

0.03799 

0.00877 

0.00023 

0.00877 

0.00023 

0.37000 

0.03799 

0.00816 

0.00067 

0.00816 

JM AiMUIT 
V • KWV O f 

0 . 38000 

0.03799 

0.00793 

0.00036 

0.00793 

0.00036 

0.39000 

0.03799 

0.00729 

-0.00039 

0.00729 

-0.00039 

0 . 40000 

0.03799 

0.00614 

-0.00073 

0.00614 

-0.00073 

0.41000 

0.03799 

0.00614 

-0.00027 

0.00514 

-0.00027 

0.42000 

0.03799 

0.00488 

0.00049 

0.00488 

0.00049 

0.43000 

0.03799 

0.00520 

0.00070 

0.00520 

0.00070 

0.44000 

0.03799 

0.00539 

0.00016 

0.00539 

0.00016 

0.46000 

0.03799 

0.00490 

-0.00057 

0.00490 

-0.00057 

0.46000 

0.03799 

0.00400 

-0.00067 

0.00400 

-0.00067 

0.47000 

0.03799 

0.00342 

-0.00003 

0.00342 

-0.00003 

0.48000 

0.03799 

0.00358 

0.00064 

0.00358 

0.00064 

0.49000 

0.03799 

0.00412 

0.00061 

0.00412 

0.00061 

0.50000 

0.03799 

0.00426 

-0.00009 

0.00426 

-0.00009 

0.61000 

0.03799 

0.00371 

-0.00069 

0.00371 

-0.00069 

0.52000 

0.03799 

0.00293 

-0.00053 

0.00293 

-0.00053 

0.53000 

0.03799 

0.00266 

0.00021 

0.00266 

0.00021 

0.54000 

0.03799 

0.00309 

0.00072 

0.00309 

0.00072 

0.56000 

0.03799 

0.00366 

0.00044 

0.00366 

0.00044 

0.56000 

0.03799 

0.00366 

-0.00032 

0.00366 

-0.00032 

0.57000 

0.03799 

0.00303 

-0.00073 

0.00303 

-0.00073 

0.58000 

0.03799 

0.00239 

-0.00034 

0.00239 

-0.00034 
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0.59000 

0.03799 

0.00240 

0.00043 

0.00240 

0.00043 

0.60000 

0.03799 

0.00301 

0.00072 

0.00301 

0.00072 

0.61000 

0.03799 

0.00361 

0.00022 

0.00351 

0.00022 

0.62000 

0.03799 

0.00333 

- 0.00052 

0.00333 

- 0.00062 

0.63000 

0.03799 

0.00265 

- 0.00069 

0.00265 

- 0.00069 

0.64000 

0.03799 

0.00221 

- 0.00011 

0.00221 

- 0.00011 

0.65000 

0.03799 

0.00249 

0.00060 

0.00249 

0.00060 

0.66000 

0.03799 

0.00319 

0.00066 

0.00319 

0.00066 

0.67000 

0.03799 

0.00356 

- 0.00001 

0.00356 

- 0.00001 

0 . 68000 

0.03799 

0.00321 

- 0.00066 

0.00321 

- 0.00066 

0.69000 

0.03799 

0.00257 

- 0.00058 

0.00257 

- 0.00058 

0 . 70000 

0.03799 

0.00236 

0.00014 

0.00236 

0.00014 

0.71000 

0.03799 

0.00288 

0.00070 

0.00288 

0.00070 

0.72000 

0.03799 

0.00360 

0.00050 

0.00360 

0.00060 

0.73000 

0.03799 

0.00381 

- 0.00025 

0.00381 

- 0.00025 

0.74000 

0.03799 

0.00336 

- 0.00073 

0.00335 

- 0.00073 

0.76000 

0.03799 

0.00282 

- 0.00040 

0.00282 

- 0.00040 

0.76000 

0.03799 

0.00290 

0.00036 

0.00290 

0.00036 

0.77000 

0.03799 

0.00362 

0.00073 

0.00362 

0.00073 

0 . 78000 . 

0.03799 

0.00432 

0.00030 

0.00432 

0.00030 

0.79000 

0.03799 

0.00438 

- 0.00046 

0.00438 

- 0.00046 

0 . 80000 

0.03799 

0.00390 

- 0.00071 

0.00390 

- 0.00071 

0.81000 

0.03799 

0.00359 

- 0.00018 

0.00369 

- 0.00018 

0.82000 

0.03799 

0.00401 

0.00055 

0.00401 

0.00066 

0.83000 

0.03799 

0.00494 

0.00068 

0.00494 

0.00068 

0.84000 

0.03799 

0.00562 

0.00006 

0.00562 

0.00006 

0.85000 

0.03799 

0.00563 

- 0.00062 

0.00663 

- 0.00062 

0.86000 

0.03799 

0.00532 

- 0.00062 

0.00632 

- 0.00062 

0.87000 

0.03799 

0.00544 

0.00006 

0.00644 

0.00006 

0.88000 

0.03799 

0.00635 

0.00067 

0.00636 

0.00067 

0 . 89000 

0.03799 

0.00763 

0.00064 

0.00763 

0.00064 

0 . 90000 

0.03799 

0.00866 

- 0.00017 

0.00865 

- 0.00017 

0.91000 

0.03799 

0.00902 

- 0.00062 

0.00902 

- 0.00062 

0.92000 

0.03799 

0.00979 

- 0.00032 

0.00979 

- 0.00032 

0.93000 

0.03799 

0.01186 

0.00005 

0.01186 

0.00006 

0.94000 

0.03799 

0.01523 

- 0.00113 

0.01624 

- 0.00109 

0.95000 

0.03799 

0.01836 

- 0.00485 

0.01838 

- 0.00474 

0 . 96000 

0.03799 

0.01948 

- 0.00990 

0.01966 

- 0.00973 

0.97000 

0.03799 

0.01872 

- 0.01382 

0.01888 

- 0.01360 

0.98000 

0.03799 

0.01813 

- 0.01561 

0.01836 

- 0.01536 

1.26620 

0.00168 

- 0.30135 

- 0.21662 

- 0.22326 

- 0.29566 
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SOURCE DISTRIBUTION REPRESENTING CIRCULAR BODY 
NSOR =71 NPRT = 2 



INPUT BODY GEOMETRY AT XFC( 1) 


7 . 19000 


XRC (Y) 

YRC (Z) 

0.0000 

-1 .6675 

0.26S8 

-1 .6433 

0.5437 

-1.5708 

0.8338 

-1.4621 

1.6433 

-0.7613 

1.9575 

-0.2973 

2.9000 

0 . 0000 

1.8850 

0.3142 

0.8458 

1.2083 

0.5437 

1.5104 

0.2900 

1.7279 

0 . 0000 

1.8367 


BODY AXIS SHIFT, ZZ = 0.0313693 
CIRCLE RADIUS, RZ = 1.9581980 
MAPPED AND INPUT CORNER COORDINATES 

MAPPED . . . INPUT . . . 


I 

YRT 

ZRT 

YR 

ZR 

1 

0 . 00000 

-1.68264 

0.00000 

-1.66750 

2 

0.26594 

-1.65941 

0.26883 

-1.64333 

3 

0.64412 

-1.58802 

0.64376 

-1.67083 

4 

0.83281 

-1.47680 

0.83376 

-1.46208 

5 

1.64575 

-0.77882 

1.64333 

-0.76125 

6 

1.98244 

-0.34142 

1.96760 

-0.29725 

7 

2.90258 

-0.01739 

2.90000 

0.00000 

8 

1.88626 

0.31646 

1.88500 

0.31417 

9 

1.03457 

1.03081 

0.84683 

1.20833 

10 

0.54531 

1.51249 

0.64375 

1.61042 

11 

0.29004 

1.72881 

0.29000 

1.72792 

12 

0.00000 

1.83667 

0.00000 

1.83667 


(j) Page 1 0 


Figure 26.- Continued. 



ITERATION NO. 


12 ***** 


SERIES COEFFICIENTS 
N AN 


1 

- 0.2646261 

2 

- 0.0172847 

3 

0.1918155 

4 

- 0.0350858 

5 

- 0.0334423 

6 

- 0.0042402 

7 

- 0.0029773 

8 

0.0068138 

9 

0.0072526 

10 

- 0.0043737 

11 

- 0.0047305 

12 

- 0.0004754 

13 

- 0.0003627 

14 

0.0052962 

15 

0.0062696 

16 

- 0.0021943 

17 

- 0.0025676 

18 

- 0.0022914 

19 

- 0.0016160 

20 

0.0021578 

21 

0.0022285 

22 

- 0.0013290 

23 

- 0.0002803 

24 

- 0.0022588 

25 

- 0.0006178 

26 

0.0014533 

27 

0.0019017 

28 

- 0.0001131 

29 

0.0003714 

30 

- 0.0014244 


(k) Page 1 1 


Figure 26.- Continued. 
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CORRESPONDING POINTS ON BODY AND CIRCLE FOR CHINE FOR X * 


I 

/ 

THETA 

- CIRCLE - 
Y 

Z 

-/ 

Y 

— BODY - 
Z 

BETA 

1 

0 . 0000 

0.0000 

-1.9582 

0.0000 

-1.6826 

0.0000 

2 

6.0000 

0.1707 

-1.9607 

0.1711 

-1.6737 

6.8381 

3 

10.0000 

0.3400 

-1.9284 

0.3421 

-1.6436 

11.7570 

4 

15.0000 

0.5068 

-1.8916 

0.6154 

-1.6971 

17.8866 

5 

20.0000 

0.6697 

-1.8401 

0.6872 

-1.6388 

24.0660 

6 

25.0000 

0.8276 

-1.7747 

0.8420 

-1.4705 

29.7952 

7 

30.0000 

0.9791 

-1.6968 

0.9784 

-1.3587 

35.7693 

8 

35.0000 

1.1232 

-1.6041 

1.1425 

-1.2137 

43.2706 

9 

40 . 0000 

1.2687 

-1.5001 

1.3124 

-1.0708 

60.7885 

10 

45 .0000 

1.3847 

-1.3847 

1.4877 

-0.9134 

58.4505 

11 

50.0000 

1.5001 

-1.2587 

1.6615 

-0.7724 

64.9362 

12 

55.0000 

1.6041 

-1.1232 

1.7783 

-0.5693 

72.2488 

13 

60 . 0000 

1.6958 

-0.9791 

2.0147 

-0.3187 

81.0097 

14 

65 .0000 

1.7747 

-0.8276 

2.3058 

-0.1997 

85.0503 

IS 

70 . 0000 

1.8401 

-0.6697 

2.5057 

-0.1466 

86.6746 

16 

75.0000 

1.8915 

-0.5068 

2.6668 

-0.0888 

88.0928 

17 

80.0000 

1.9284 

-0.3400 

2.7864 

-0.0656 

88.8583 

18 

85.0000 

1.9607 

-0.1707 

2.8650 

-0.0272 

89.4666 

19 

90.5845 

1.9581 

0.0200 

2.9026 

-0.0174 

89.6667 

20 

95.0000 

1.9507 

0,1707 

2.8786 

-0.0097 

89.6066 

21 

100.0000 

1.9284 

0.3400 

2.8093 

0.0124 

90.2621 

22 

105.0000 

1.8915 

0.5068 

2.7014 

0.0454 

90.9636 

23 

110.0000 

1.8401 

0.6697 

2.5560 

0.0898 

92.0130 

24 

115.0000 

1.7747 

0.8276 

2.3770 

0.1475 

93.5518 

25 

120.0000 

1.6958 

0.9791 

2.1556 

0.2104 

96.5766 

26 

125.0000 

1.6041 

1.1232 

1.8721 

0.3243 

99.8284 

27 

130.0000 

1.5001 

1.2587 

1.6125 

0.6221 

107.9414 

28 

136.0000 

1.3847 

1.3847 

1.4062 

0.7102 

116.7947 

29 

140.0000 

1.2587 

1.5001 

1.1901 

0.8865 

126.6817 

30 

145.0000 

1.1232 

1.6041 

0.9895 

1.0742 

137.3611 

31 

150.0000 

0.9791 

1.6968 

0.8048 

1.2616 

147.2691 

32 

155.0000 

0.8276 

1.7747 

0.6430 

1.4196 

166.6313 

33 

160.0000 

0.6697 

1.8401 

0.4953 

1.6666 

162.3392 

34 

165.0000 

0.5068 

1.8916 

0.3587 

1,6764 

167.9216 

35 

170.0000 

0.3400 

1.9284 

0.2422 

1.7517 

172.1283 

36 

175.0000 

0.1707 

1.9607 

0.1056 

1.8008 

176.6437 

37 

180.0000 

0.0000 

1.9682 

0 . 0000 

1.8367 

180.0000 

38 

185.0000 

-0.1707 

1.9507 

-0.1056 

1 . 8008 

183.3663 

39 

190 . 0000 

-0.3400 

1.9284 

-0.2422 

1.7617 

187.8717 

40 

195.0000 

-0.5068 

1.8916 

-0.3587 

1.6764 

192.0785 

41 

200.0000 

-0.6697 

1.8401 

-0.4953 

1.5656 

197.6608 

42 

205.0000 

-0.8276 

1.7747 

-0.6430 

1.4195 

204.3687 

43 

210.0000 

-0.9791 

1.6958 

-0.8048 

1.2516 

212.7410 

44 

215.0000 

-1.1232 

1.6041 

-0.9895 

1.0742 

222.6489 

45 

220.0000 

-1.2587 

1.6001 

-1.1901 

0.8865 

233.3184 

46 

225.0000 

-1.3847 

1.3847 

-1.4062 

0.7102 

243.2053 

47 

230.0000 

-1.6001 

1.2687 

-1 .6125 

0.6221 

262.0686 

48 

235.0000 

-1.6041 

1.1232 

-1.8721 

0.3243 

260.1716 

49 

240.0000 

-1.6958 

0.9791 

-2.1566 

0.2104 

264.4246 

50 

245.0000 

-1.7747 

0.8276 

-2.3770 

0.1476 

266.4482 

51 

250.0000 

-1.8401 

0.6697 

-2.6660 

0.0898 

267.9870 

52 

265.0000 

-1.8915 

0.5068 

-2.7014 

0.0454 

269.0366 

53 

260.0000 

-1.9284 

0.3400 

-2.8093 

0.0124 

269.7479 

64 

265.0000 

-1.9607 

0.1707 

-2.8786 

-0.0097 

270.1934 

55 

269.4155 

-1.9581 

0.0200 

-2.9026 

-0.0174 

270.3433 

66 

276.0000 

-1.9507 

-0.1707 

-2.8650 

-0.0272 

270.6444 

67 

280.0000 

-1.9284 

-0.3400 

-2.7864 

-0.0666 

271.1417 

58 

285.0000 

-1.8915 

-0.5068 

-2.6658 

-0.0868 

271.9072 

69 

290.0000 

-1.8401 

-0.6697 

-2.6057 

-0.1456 

273.3254 

60 

295.0000 

-1.7747 

-0.8276 

-2.3058 

-0.1997 

274.9497 

61 

300.0000 

-1.6958 

-0.9791 

-2.0147 

-0.3187 

278.9904 


7.19000 

/ 


* CHINE SEPARATION POINT 


• CHINE SEPARATION POINT 


(I) Page 12 


Figure 26.- Continued. 



62 

305.0000 

-1.6041 

-1.1232 

63 

310.0000 

-1.5001 

-1.2587 

64 

315.0000 

-1.3847 

-1.3847 

65 

320.0000 

-1.2587 

-1.5001 

66 

325.0000 

-1.1232 

-1.6041 

67 

330.0000 

-0.9791 

-1.6958 

68 

335.0000 

-0.8276 

-1.7747 

69 

340.0000 

-0.6697 

-1.8401 

70 

346.0000 

-0.5068 

-1.8915 

71 

3S0.0000 

-0.3400 

-1.9284 

72 

355.0000 

-0.1707 

-1.9507 

73 

360.0000 

0 . 0000 

-1.9682 


4 ^ 


-1.7783 
-1.6S16 
-1.4877 
-1.3124 
-1.1426 
-0.9784 
-0.8420 
-0.6872 
-0 .5154 
-0.3421 
-0.1711 
0.0000 


-0.B693 

-0.7724 

-0.9134 

-1.0708 

-1.2137 

-1.3B87 

-1.4705 

-1.6388 

-1.B971 

-1.6436 

-1.6737 

-1.6826 


287.7513 

296.0648 

301.6496 

309.2116 

316.7296 

324.2408 

330.2048 

335.9340 

342.1136 

348.2430 

354.1619 

0.0000 


(m) Page 13 

Figure 26.- Continued. 



INPUT BODY GEOMETRY AT XFC( 2) 


13.66000 


XRC (Y) 

YRC (Z) 

0.0000 

-2.3200 

0.4108 

-2.2958 

0.7854 

-2.1992 

1.2083 

-2.0300 

2.3563 

-1.1358 

2.9000 

-0.4350 

4.1325 

0.0000 

2.8758 

0.4471 

1.1479 

2.0300 

0.7975 

2.2958 

0.3988 

2.5133 

0 . 0000 

2.5979 


BODY AXIS SHIFT, ZZ =-0.0303943 
CIRCLE RADIUS, RZ = 2.8173276 
MAPPED AND INPUT CORNER COORDINATES 

MAPPED... INPUT... 


I YRT 


ZRT 


YR ZR 


1 

0.00000 

-2.31520 

2 

0.41097 

-2.28988 

3 

0.78591 

-2.19382 

4 

1.20788 

-2.02176 

5 

2.35721 

-1.12858 

6 

2.93405 

-0.45531 

7 

4.13458 

0.00735 

8 

2.89800 

0.47194 

9 

1.15134 

2.03761 

10 

0.79788 

2.29922 

11 

0.39795 

2.51455 

12 

0 . 00000 

2.59792 


0.00000 -2.32000 
0.41083 -2.29583 
0.78542 -2.19917 
1.20833 -2.03000 
2.35625 -1 . 13583 
2.90000 -0.43500 
4.13250 0.00000 
2.87583 0.44708 
1.14792 2.03000 
0.79750 2.29583 
0.39875 2.61333 
0.00000 2.59792 


(n) Page 14 


Figure 26.- Continued. 



•**•* ITERATION NO. 


13 


SERIES COEFFICIENTS 
N AN 


1 

- 0.2611494 

2 

- 0. 0208124 

3 

0.1685687 

4 

- 0.0170018 

5 

- 0.0405029 

8 

- 0.0006268 

7 

0.0027656 

8 

0.0026457 

9 

0.0061041 

10 

- 0.0036377 

11 

- 0.0066124 

12 

- 0.0004192 

13 

0.0009768 

14 

0.0026675 

15 

0.0055238 

16 

- 0.0007471 

17 

- 0.0041913 

18 

- 0.0011172 

19 

- 0.0006133 

20 

0.0005546 

21 

0.0019569 

22 

- 0.0005039 

23 

- 0.0016407 

24 

- 0.0005206 

25 

- 0.0003263 

26 

0.0003213 

27 

0.0017617 

28 

- 0.0001255 

29 

- 0.0006746 

30 

- 0.0000636 


(o) Page 15 


Figure 26.- Continued. 




CORRESPONDING POINTS 

ON BODY , 

AND CIRCLE 


/ 

-- CIRCLE - 


— / 

I 

1IILTA 

Y 

Z 

Y 

1 

0.0000 

0.0000 

-2.8173 

0 .0000 

2 

6.0000 

0.2455 

-2.8066 

0.2671 

3 

10.0000 

0.4892 

-2.7745 

0.5041 

4 

15.0000 

0.7292 

-2.7213 

0.7498 

6 

20.0000 

0.9636 

-2.6474 

0.9977 

6 

25.0000 

1.1907 

-2.6634 

1.2267 

7 

30.0000 

1.4087 

-2.4399 

1.4430 

8 

35.0000 

1.6160 

-2.3078 

1.6883 

9 

40.0000 

1.8109 

-2.1682 

1.9444 

10 

45.0000 

1.9922 

-1.9922 

2.1966 

11 

60.0000 

2.1582 

-1.8109 

2.4222 

12 

65.0000 

2.3078 

-1.6160 

2.6145 

13 

60.0000 

2.4399 

-1.4087 

2.9213 

14 

65.0000 

2.6534 

-1.1907 

3.3070 

16 

70.0000 

2.6474 

-0.9636 

3.6922 

16 

76.0000 

2.7213 

-0.7292 

3.8111 

17 

80.0000 

2.7745 

-0.4892 

3.9818 

18 

85.0000 

2.8066 

-0.2465 

4.0882 

19 

89.9071 

2.8173 

-0.0046 

4.1346 

20 

95.0000 

2.8066 

0.2455 

4.0860 

21 

100.0000 

2.7746 

0.4892 

3.9766 

22 

106.0000 

2.7213 

0.7292 

3.8092 

23 

110.0000 

2.6474 

0.9636 

3.6889 

24 

116.0000 

2.5634 

1.1907 

3.3197 

25 

120.0000 

2.4399 

1.4087 

2.9574 

26 

125.0000 

2.3078 

1.6180 

2,6050 

27 

130.0000 

2.1582 

1.8109 

2.3241 

28 

135.0000 

1.9922 

1.9922 

2.0461 

29 

140.0000 

1.8109 

2.1582 

1.7761 

30 

145.0000 

1.6160 

2.3078 

1.6192 

31 

150.0000 

1.4087 

2.4399 

1.2734 

32 

155.0000 

1.1907 

2.5634 

1.0488 

33 

160.0000 

0.9636 

2.6474 

0.8248 

34 

165.0000 

0.7292 

2.7213 

0.6084 

35 

170.0000 

0.4892 

2.7746 

0.4025 

36 

175.0000 

0.2455 

2.8066 

0.1964 

37 

180.0000 

0.0000 

2.8173 

0 . 0000 

38 

185.0000 

-0.2466 

2 . 8066 

-0.1964 

39 

190.0000 

-0.4892 

2.7746 

-0.4025 

40 

196.0000 

-0.7292 

2.7213 

-0.6084 

41 

200.0000 

-0.9636 

2.6474 

-0.8248 

42 

205.0000 

-1.1907 

2.5534 

-1.0488 

43 

210.0000 

-1.4087 

2.4399 

-1.2734 

44 

216.0000 

-1.6160 

2.3078 

-1.5192 

45 

220.0000 

-1.8109 

2.1682 

-1.7761 

46 

226.0000 

-1.9922 

1.9922 

-2.0461 

47 

230.0000 

-2.1682 

1.8109 

-2.3241 

48 

235.0000 

-2.3078 

1.6160 

-2.6050 

49 

240.0000 

-2.4399 

1.4087 

-2.9674 

60 

246.0000 

-2.5534 

1.1907 

-3.3197 

61 

250.0000 

-2.6474 

0.9636 

-3.6889 

52 

255.0000 

-2.7213 

0.7292 

-3.8092 

53 

260.0000 

-2.7745 

0.4892 

-3.9766 

64 

265.0000 

-2.8066 

0.2465 

-4.0860 

56 

270.0929 

-2.8173 

-0.0046 

-4 . 1346 

56 

276.0000 

-2.8066 

-0.2455 

-4.0882 

67 

280.0000 

-2.7745 

-0.4892 

-3.9818 

58 

285.0000 

-2.7213 

-0.7292 

-3.8111 

69 

290.0000 

-2.6474 

-0.9636 

-3.5922 

60 

295.0000 

-2.5534 

-1.1907 

-3.3070 

61 

300.0000 

-2.4399 

-1.4087 

-2.9213 


FOR CHINE FOR X = 13.66000 


-- BODY / 

Z BETA 

-2.3162 0.0000 

-2.3071 6.3698 

-2.2729 12.6080 

-2.2067 18.7736 

-2.1180 26.2218 

-2.0108 31.3864 

-1.8614 37.9337 

-1.6667 46.6698 

-1.4610 63.0782 

-1.2609 60.1432 

-1.0637 66.2923 

-0.7926 73.1366 

-0.4646 80.9648 

-0.2760 86.2287 

-0.1894 86.9814 

-0.1037 88.4416 

-0.0492 89.2918 

aaQ 1 on oqai 

0^0073 90 ! 1018 * CHINE SEPARATION POINT 

0.0215 90.3016 

0.0663 90.9410 

0.1194 91.7949 

0.2039 93.2613 

0.2928 96.0402 

0.4376 98.4164 


0.6992 105.0266 

0.9706 112.6667 

1.2182 120.7688 

1.4690 129.6930 

1.7040 138.2816 

1.9278 146.6636 

2.1223 163.7031 

2.2816 160.1248 

2.4168 166.8712 

2.6129 170.8990 

2.6734 176.6360 

2.6979 180.0000 

2.5734 184.3641 

2.6129 189.1010 

2.4168 194.1286 

2.2816 199.8763 

2.1223 206.2969 

1.9278 213.4464 

1.7040 221.7186 

1.4690 230.4070 

1.2182 239.2312 

0.9706 247.3333 

0.6992 264.9746 

0.4376 261.6836 

0.2928 264.9698 

0.2039 266.7487 

0.1194 268.2061 

0.0653 269.0590 

0.0216 269.6986 

0.0073 269.8982 • CHINE SEPARATION POINT 

-0.0081 270.1139 

-0.0492 270.7082 

-0.1037 271.6683 

-0.1894 273.0186 

-0.2760 274.7713 

-0.4645 279.0362 
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62 

305.0000 

- 2.3078 

- 1.6160 

- 2.6145 

- 0.7925 

286.8636 

63 

310.0000 

- 2.1582 

- 1.8109 

- 2.4222 

- 1.0637 

293.7677 

64 

315.0000 

- 1.9922 

- 1.9922 

- 2.1966 

- 1.2609 

299.8669 

65 

320.0000 

- 1.8109 

- 2.1582 

- 1.9444 

- 1.4610 

306.9219 

66 

325.0000 

- 1.6160 

- 2.3078 

- 1.6883 

- 1.6557 

314.4402 

67 

330.0000 

- 1.4087 

- 2.4399 

- 1.4430 

- 1.8514 

322.0663 

68 

335.0000 

- 1.1907 

- 2.5534 

- 1.2267 

- 2.0108 

328.6136 

69 

340.0000 

- 0.9636 

- 2.6474 

- 0.9977 

- 2.1180 

334.7783 

70 

345.0000 

- 0.7292 

- 2.7213 

- 0.7498 

- 2.2067 

341.2264 

71 

350.0000 

- 0.4892 

- 2.7745 

- 0.5041 

- 2.2729 

347.4960 

72 

355.0000 

- 0.2455 

- 2.8066 

- 0.2671 

- 2.3071 

363.6402 

73 

360.0000 

0 . 0000 

- 2.8173 

0.0000 

- 2.3162 

0.0000 
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INPUT BODY GEOMETRY AT XFC ( 3) = 19.94000 


XRC (Y) 

YRC (Z) 

0 . 0000 

-2.6524 

0.4606 

-2.6121 

0.9072 

-2.4925 

1.3262 

-2.2971 

2.4111 

-1.2307 

3.1590 

-0.4827 

4.3500 

0.0000 

3.1590 

0.4827 

1.3262 

2.2971 

0.9072 

2.4925 

0.4606 

2.6121 

0.0000 

2.6524 


BODY AXIS SHIFT, 22 = 0.0002265 
CIRCLE RADIUS, R2 = 3.0342150 
MAPPED AND INPUT CORNER COORDINATES 

MAPPED. . . INPUT. . . 


I 

YRT 

2RT 

YR 

ZR 

1 

0 . 00000 

-2.68137 

0.00000 

-2.65240 

2 

0.46035 

-2.61169 

0.46060 

-2.61210 

3 

0.90681 

-2.49268 

0.90719 

-2.49247 

4 

1.32411 

-2.29426 

1.32622 

-2.29709 

5 

2.41363 

-1.23347 

2.41107 

-1.23073 

6 

3.17855 

-0.50470 

3.15905 

-0.48274 

7 

4.35134 

-0.00256 

4.35000 

0.00000 

8 

3.19524 

0.49184 

3.15905 

0.48274 

9 

1.32518 

2.29484 

1.32622 

2.29709 

10 

0.90717 

2.49336 

0.90719 

2.49247 

11 

0.45997 

2.61266 

0.46060 

2.61210 

12 

0.00000 

2.65240 

0.00000 

2.66240 
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ITERATION NO. 


12 


SERIES COEFFICIENTS 
N AN 


1 

- 0 . 2313311 

2 

- 0 . 0004480 

3 

0 . 1544752 

4 

0.0002782 

B 

- 0.0517554 

6 

- 0.0000627 

7 

- 0 . 0007065 

8 

0.0000508 

9 

0.0037163 

10 

0.0000468 

11 

- 0.0030484 

12 

- 0.0001696 

13 

0.0028445 

14 

0.0000606 

15 

0.0031767 

16 

0.0001042 

17 

- 0.0038550 

18 

- 0.0001290 

19 

0.0005987 

20 

0.0000396 

21 

0.0001602 

22 

0.0000740 

23 

- 0.0011882 

24 

- 0.0001046 

25 

0.0012832 

26 

0.0000243 

27 

0.0007353 

28 

0.0000586 

29 

- 0.0011114 

30 

- 0.0000674 
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Figure 26.- Continued. 




CORRESPONDING POINTS 

ON BODY 

AND CIRCLE 

FOR CHINE 

FOR X = 19.94000 




/ 

— CIRCLE - 


-/ 

BODY - 




I 

THETA 

Y 

Z 

Y 

2 

BETA 



1 

0.0000 

0.0000 

-3.0342 

0 . 0000 

-2.6514 

0 . 0000 



2 

5 . 0000 

0.2644 

-3.0227 

0.2643 

-2.6397 

6.7167 



3 

10.0000 

0.5269 

-2.9881 

0.6253 

-2.6990 

11.4254 



4 

15.0000 

0.7853 

-2.9308 

0.7835 

-2.5344 

17.1781 



S 

20.0000 

1.0378 

-2.8512 

1.0327 

-2.4424 

22.9185 



6 

25.0000 

1.2823 

-2.7499 

1.2689 

-2.3288 

28.5861 



7 

30.0000 

1.6171 

-2.6277 

1.4728 

-2.1679 

34.1909 



8 

35.0000 

1.7404 

-2.4866 

1.6961 

-1.9364 

41.2152 



9 

40.0000 

1.9504 

-2.3243 

1.9525 

-1.6937 

49.0602 



10 

45.0000 

2.1455 

-2.1455 

2.2090 

-1.4334 

57.0207 



11 

60.0000 

2.3243 

-1.9604 

2.4875 

-1.1645 

64.9142 



12 

65.0000 

2.4855 

-1.7404 

2.7654 

-0.8922 

72.0686 



13 

60.0000 

2.6277 

-1.5171 

3.0670 

-0.6856 

79.1904 



14 

65.0000 

2.7499 

-1.2823 

3.4464 

-0.3642 

83.9663 



15 

70.0000 

2.8512 

-1.0378 

3.7672 

-0.2466 

86.2447 



16 

75.0000 

2.9308 

-0.7853 

3.9936 

-0.1475 

87.8854 



17 

80.0000 

2.9881 

-0.6269 

4.1786 

-0.0735 

88.9921 



18 

85.0000 

3.0227 

-0.2644 

4.3017 

-0.0231 

89.6927 



19 

89.9833 

3.0342 

-0.0009 

4.3513 

-0.0026 

89.9662 * CHINE 

SEPARATION 

POINT 

20 

95.0000 

3.0227 

0.2644 

4.3007 

0.0186 

90.2479 



21 

100.0000 

2.9881 

0.5269 

4.1769 

0.0688 

90.9433 



22 

105.0000 

2.9308 

0.7863 

3.9908 

0.1435 

92.0699 



Z3 

110.0000 

2.8612 

1.03/8 

3 . 7538 

0.2423 

93.6925 



24 

115.0000 

2.7499 

1.2823 

3.4381 

0.3620 

96.0100 



25 

120.0000 

2.6277 

1.5171 

3.0626 

0.6893 

100.8908 



26 

125.0000 

2.4855 

1.7404 

2.7669 

0.8928 

107.9508 



27 

130.0000 

2.3243 

1.9604 

2.4860 

1.1628 

116.0763 



28 

135.0000 

■2.1456 

2.1455 

2.2059 

1.4339 

123.0261 



29 

140.0000 

1.9504 

2.3243 

1.9498 

1.6946 

130.9949 



30 

145.0000 

1.7404 

2.4856 

1 .6939 

1.9385 

138.8526 



31 

150.0000 

1.6171 

2.6277 

1.4721 

2.1699 

146.8462 



32 

155.0000 

1.2823 

2.7499 

1.2686 

2.3304 

161.4374 



33 

160.0000 

1.0378 

2.8612 

1.0327 

2.4436 

157.0900 



34 

165.0000 

0.7853 

2.9308 

0.7832 

2.6362 

162.8333 



36 

170.0000 

0.5269 

2.9881 

0.6253 

2.6000 

168.5786 



36 

175.0000 

0.2644 

3.0227 

0.2641 

2.6404 

174.2877 



37 

180.0000 

0.0000 

3.0342 

0.0000 

2.6524 

180.0000 



38 

185.0000 

-0.2644 

3.0227 

-0.2641 

2.6404 

185.7123 



39 

190.0000 

-0.6269 

2.9881 

-0.5253 

2.6000 

191.4214 



40 

195.0000 

-0.7853 

2.9308 

-0.7832 

2.5352 

197.1667 



41 

200.0000 

-1.0378 

2.8612 

-1.0327 

2.4436 

202.9101 



42 

205.00CC 

-1.2823 

2.7439 

-1.2636 

2.3304 

208.6626 



43 

210.0000 

-1.5171 

2.6277 

-1.4721 

2.1699 

214.1639 



44 

215.0000 

-1.7404 

2.4865 

-1.6939 

1.9385 

221.1475 



45 

220.0000 

-1.9504 

2.3243 

-1.9498 

1.6946 

229.0062 



46 

225.0000 

-2.1455 

2.1455 

-2.2059 

1.4339 

236.9749 



47 

230.0000 

-2.3243 

1.9504 

-2.4850 

1.1628 

244.9237 



48 

235.0000 

-2.4855 

1.7404 

-2.7559 

0.8928 

252.0493 



49 

240.0000 

-2.6277 

1.6171 

-3.0626 

0.6893 

269.1093 



50 

246.0000 

-2.7499 

1.2823 

-3.4381 

0.3620 

263.9901 



61 

250.0000 

-2.8512 

1.0378 

-3.7538 

0.2423 

266.3075 



62 

255.0000 

-2.9308 

0.7853 

-3.9908 

0.1435 

267.9401 



63 

260.0000 

-2.9881 

0.6269 

-4.1769 

0.0688 

269.0667 



64 

265.0000 

-3.0227 

0.2644 

-4.3007 

0.0186 

289.7521 



65 

270.0167 

-3.0342 

-0.0009 

-4.3513 

-0.0026 

270.0338 * CHINE 

SEPARATION 

POINT 

66 

275.0000 

-3.0227 

-0.2646 

-4.3017 

-0.0231 

270.3073 



67 

280.0000 

-2.9881 

-0.5269 

-4.1786 

-0.0735 

271.0079 



68 

285.0000 

-2.9308 

-0.7853 

-3.9936 

-0.1475 

272.1147 



69 

290.0000 

-2.8612 

-1.0378 

-3.7572 

-0.2466 

273.7663 



60 

295.0000 

-2.7499 

-1.2823 

-3.4454 

-0.3642 

276.0337 



61 

300.0000 

-2.6277 

-1.5171 

-3.0670 

-0.5856 

280.8096 
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62 

305.0000 

- 2.4855 

- 1 .7404 

- 2.7554 

- 0.8922 

287.9416 

63 

310.0000 

- 2.3243 

- 1.9504 

- 2.4875 

- 1.1645 

295.0858 

64 

315.0000 

- 2.1455 

- 2.1455 

- 2.2090 

- 1.4334 

302.9794 

65 

320.0000 

- 1.9504 

- 2.3243 

- 1.9525 

- 1.6937 

310.9398 

66 

325.0000 

- 1.7404 

- 2.4855 

- 1.6961 

- 1.9364 

318.7849 

67 

330.0000 

- 1.5171 

- 2.6277 

- 1.4728 

- 2.1679 

326.8091 

68 

335.0000 

- 1.2823 

- 2.7499 

- 1.2689 

- 2.3288 

331.4139 

69 

340.0000 

- 1.0378 

- 2.8512 

- 1.0327 

- 2.4424 

337.0815 

70 

345.0000 

- 0.7853 

- 2.9308 

- 0.7835 

- 2.5344 

342.8219 

71 

350.0000 

- 0.5269 

- 2.9881 

- 0.6253 

- 2.5990 

348.5746 

72 

355.0000 

- 0.2644 

- 3.0227 

- 0.2643 

- 2.6397 

354.2834 

73 

360.0000 

0.0000 

- 3.0342 

0 » 0000 

- 2.6514 

0 . 0000 
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CALCULATED COEFFICIENTS FOR NUMERICAL MAPPING, MNFC = 30 MXFC = 3 

XFC 

7.190 


BODY AXIS SHIFT, ZZ = 0.0313693 
CIRCLE RADIUS, RZ = 1.9681980 


XRC (Y) 



0 .000008+00 
2 . 902588+00 

2.66938E-01 
1 . 88626E+00 

6.44120E-01 
1 .03467E+00 

8 . 32810E-01 
6.45313E-01 

1 . 64676E+00 
2 . 90037E-01 

1 . 98244E+00 
-8 . 83767E-08 

YRC 

(Y) 

-1 .682048*00 
-1 .73853E-02 

>1 .65941E+00 
3 . 16462E-01 

-1 .58802E+00 
1 .030818+00 

-1 . 47680E+00 
1 . 61249E+00 

-7 . 78818E-01 
1.72881E+00 

-3 . 41417E-01 
1 . 83667E+00 

THC 

(RAD.) 

0 . 00000E+00 
1 .58119E+00 

1 .35886E-01 
2 . 17743E+00 

2 . 76273E-01 
2.610146+00 

4 . 30621E-01 
2.76318E+00 

8 . 68877E-01 
2 . 93224E+00 

1 . 03771E+00 
3 . 141S9E+00 

AFC 

-2 .846268-01 
-2 . 97731E-03 
-3.62738E-04 
-1 .61S99E-03 
-6 . 17825E-04 

-1 .72847E-02 
6 . 81383E-03 
6 . 29017E-03 
2 . 16777E-03 
1 . 45334E-03 

1 .91816E-01 
7 . 26260E-03 
6 . 26960E-03 
2 . 22853E-03 
1 . 90170E-03 

-3 . 60868E-02 
-4 . 37366E-03 
-2.19430E-03 
-1 .32900E-03 
-1 . 13097E-04 

-3 . 34423E-02 
-4 . 73047E-03 
-2 . 66767E-03 
-2 . 80266E-04 
3.71434E-04 

-4 . 24022E-03 
-4 . 7S440E-04 
-2 . 29141E-03 
-2 . 26879E-03 
-1 . 42439E-03 

XFC 

13.660 






BODY AXIS 

SHIFT, ZZ =-0, 

,0303943 





CIRCLE RADIUS, RZ =2. 

,8173276 





XRC 

(Y) 

0 . 00000E+00 
A. I3468E+00 

4 . 10966E-01 
2 . 83800E+0B 

7 . 86907E-01 
1. 151348+00 

1 . 20788E+00 
7 . 97879E-01 

2 . 36721E+00 
3 . 97962E-01 

2 . 93406E+00 
-1 . 92761E-07 

YRC 

(Y) 

-2.316208+00 
7 . 34697E-03 

-2.289888+00 
4 . 71940E-01 

-2 . 19382E+00 
2 . 03761E+00 

-2.02176E+00 
2 . 29922E+00 

-1 . 12868E+00 
2 . 61466E+00 

-4 . 66311E-01 
2.69792E+00 

THC 

(RAD.) 

0 . 00000E+00 
1 . 6691 1E+00 

1.41308E-01 
2 . 107968+00 

2.74484E-01 

2.66487E+00 

4 . 28731E-01 
2 . 80304E+00 

8.44968E-01 
2 . 96900E+00 

1 . 06009E+00 
3 . 14169E+00 

AFC 

-2.61 1 49E-01 
2 , 76667E-03 
9 . 76762E-04 
“6 . 13266E-04 
- 3 . 26334E-04 

-2.08124E-02 
2 . 64 670E-03 
2.66747E-03 
6.64568E-04 
3 . 21292E-04 

1 .686698-01 
6.104 10E-03 
6 . 62376E-03 
1 .956918-03 
1 . 76171E-03 

-1 . 70018E-02 
-3.63772E-03 
-7 . 47092E-04 
-6.03931E-04 
-1 . 26529E-04 

-4.06029E-02 
-6.61244E-03 
-4 . 191278-03 
-1 . 64074E-03 
-0 . 746648-04 

-0 . 2684 IE-04 
-4.192198-04 
-1 . 1 1723E-03 
-6.20813E-04 
-6 . 36085E-06 
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XFC 

19.940 


BODY AXIS SHIFT, ZZ = 0.000226S 
CIRCLE RADIUS , RZ = 3.0342150 


XRC (Y) 

0 . 00000E+00 4 . 60352E-01 9.06812E-01 1.324UE+00 2.41363E+00 3.17855E+00 

4 . 35134E+00 3 . 19524E+00 1.32518E+00 9.07167E-01 4.59968E-01 -2.66463E-07 

YRC (Y) 

“2 . 65137E+00 -2 . 61169E+00 -2.49268E+00 -2.29426E+00 -1.23347E+00 -6.04703E-01 

-2.56398E-03 4.91845E-01 2.29484E+00 2.49336E+00 2.61266E+00 2.6S240E+00 


THC (RAD.) 

0 . 00000E+00 1 . 62870E-01 3.047S3E-01 4.58858E-01 8.49708E-01 1.07326E+00 

1 . 67051E+00 2.06304E+00 2.68209E+00 2.83668E+00 2.98887E+00 3.14159E+00 


AFC 

-2.31331E-01 -4 . 48035E-04 1.54475E-01 

-7 . 06547E-04 S.07649E-05 3.71627E-03 

2 . 84450E-03 6 . 05765E-05 3 . 17672E-03 

5 . 98718E-04 3 . 96489E-05 1 . 60234E-04 

1 .28322E-03 2.42639E-05 7.35292E-04 


2 . 78203E-04 -5.17554E-02 -6.27261E-05 
4.68376E-06 -3.04841E-03 -1.69560E-04 
1 . 04184E-04 -3 . 86505E-03 -1.29021E-04 
7.40022E-05 -1 . 18822E-03 -1.04638E-04 
6 . 86612E-06 -1.U143E-03 -6.73980E-0S 


(w) Page 23 


Figure 26.- Continued. 



-ssi- 


PARAMETERS FOR CHINE CROSS SECTION BODIES 

VORTEX LATTICE APPROX. AT FIRST 8 STATIONS TO START THE CHINE SOLUTION 

VORTEX LATTICE BOUNDARY CONDITION APPLIED AT 0.26000 OF THE EXPOSED STRAKE SPAN 
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X R DR/DX X/L RETR 

6.6000 0.1491 0.3031 0.0068 0.197E+06 


BODY SURFACE PRESSURE DISTRIBUTION 
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CP 
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1 
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2 
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7 
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8 
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0.756 

-0 . 063 

9 
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10 

0.1133 
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0 . 0000 
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11 
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0 . 0000 
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0 . 858 
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12 
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72.249 
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0.0000 
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-0 . 033 

13 
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-0.0243 

81.010 

0.7951 
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0.3934 
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14 
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0.8687 
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15 

0.1908 

-0.0111 
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0.8800 

0.6610 

-0.0056 

0.6610 
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0.0000 
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0.6222 

0.609 

0.728 

16 

0.2030 

-0.0068 

88.093 

0.9307 

0.7822 

-0.0473 

0.7837 
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0.0000 

0.4021 

0.5763 

0.682 

0.951 

17 

0.2121 

-0.0042 

88.858 

0.9692 
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18 

0.2181 
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20 
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SUMMARY OF VORTEX FIELD AT X = 1.000 
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Figure 26.- Continued. 
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SUMMARY OF VORTEX FIELD AT X = 4.600 
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SUMMARY OF VORTEX FIELD AT X = 6.00C 


H = 0.60000 
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1.13254 
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Figure 26.- Continued. 
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Figure 26.- Concluded. 
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